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ABSTRACT 


The  behavior  of  the  fully  rough  turbulent  boundary  layer  subjected 
to  favorable  pressure  gradients  was  investigated  experimentally  using  a 
porous  test  surface  composed  of  densely  packed  spheres  of  uniform  size. 
Measurements  of  profiles  of  mean  velocity,  mean  temperature  and  the  com- 
ponents of  the  Reynolds  stress  tensor  are  reported  for  both  unblown  and 
blown  layers . Stanton  numbers  were  determined  from  energy  balances  on 
the  test  surface  and  skin  friction  coefficients  from  measurements  of  the 
Reynolds  shear  stress  and  mean  velocity. 

A new  acceleration  parameter,  K , for  fully  rough  layers  is  de- 
fined and  shown  to  be  dependent  on  a characteristic  roughness  dimension 
but  independent  of  molecular  viscosity.  For  Kr  constant  and  the  blow- 
ing fraction,  F,  constant  and  greater  than  or  equal  to  zero,  it  is 
shown  that  the  fully  rough  turbulent  boundary  layer  reaches  an  equilib- 
rium state  in  which  profiles  of  the  mean  velocity  and  the  Reynolds  stress 
tensor  components  are  similar  in  the  flow  direction  and  skin  friction 
coefficient,  momentum  thickness,  boundary  layer  shape  factor,  and  the 
Clauser  shape  factor  and  pressure  gradient  parameter  all  become  constant. 
The  thermal  data  indicate  the  possibility  that  such  a layer,  with  wall 
temperature  constant,  may  approach  a state  of  equilibrium  in  the  thermal 
sense,  also.  Such  a state  would  be  characterized  by  Stanton  number  be- 
coming constant,  enthalpy  thickness  approaching  an  asynptote,  and  tem- 
perature profiles  exhibiting  similarity  in  the  flow  direction. 

For  fully  rough  turbulent  flow,  acceleration  causes  an  increase  in 
Stanton  number  compared  to  zero  pressure  gradient  values  at  the  same  en- 
thalpy thickness,  Reynolds  number,  or  position.  For  the  present  range 


of  accelerations , these  increases  were  approximately  ten  and  twenty  per- 
cent for  the  unblown  and  blown  cases , respectively.  Data  for  variable 
test  surface  temperature  cases  show  that  nondimensionally  equivalent 
positive  axial  gradients  of  freestream  velocity  and  temperature  poten- 
tial across  the  boundary  layer  have  identical  effects  on  Stanton  number. 
The  fully  rough  Stanton  number  behavior  observed  in  this  study  is  con- 
trary to  that  previously  reported  for  unblown  accelerated  smooth  wall 
layers . 

Acceleration  of  a fully  rough  layer  decreases  the  normalized  turbu- 
lent kinetic  energy  and  makes  the  turbulence  field  much  less  isotropic 
in  the  inner  region  (for  F equal  zero)  compared  to  zero  pressure  gra- 
dient fully  rough  layers.  The  values  of  the  Reynolds  shear  stress  cor- 
relation coefficients,  however,  are  unaffected  by  acceleration  or  blowing 
and  are  identical  with  values  previously  reported  for  zero  pressure  gra- 
dient smooth  and  rough  wall  flows.  Increasing  values  of  roughness 
Reynolds  number  with  acceleration  indicate  that  the  fully  rough  layer 
does  not  tend  toward  the  transitionally  rough  or  smooth  wall  state  when 
accelerated. 

An  integral  prediction  method  is  presented  which  successfully  de- 
scribes Stanton  number  behavior  in  a fully  rough  turbulent  flow  with 
variable  velocity,  wall  temperature,  and  blowing  using  only  a kernel 
function  determined  from  zero  pressure  gradient  flow  with  an  unheated 
starting  length. 
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CHAPTER  1.  INTRODUCTION 
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Although  turbulent  flow  has  been  a prime  area  of  both  theoretical  and 
experimental  research  for  the  past  several  decades,  the  present  understanding 
of  the  behavior  and  basic  mechanisms  of  tubulence  fields  is  rather  tenuous, 
at  best.  Experimental  data  on  turbulent  boundary  layers  is  either  very 
limited  or  nonexistent  for  many  classes  of  boundary  conditions.  With  the 
recent  advent  of  more  sophisticated  prediction  schemes  and  turbulence  models, 
requirements  for  more  detailed  data  on  the  turbulence  field  have  increased. 
Measurements  of  skin  friction,  Stanton  number,  and  mean  temperature  and 
velocity  fields  no  longer  provide  a sufficient  data  base  from  which  turbu- 
lence behavior  may  be  examined.  Additional  measurements  of  the  turbulence 
quantities  (fluctuations,  correlations)  are  required. 

An  experimental  study  of  the  effects  of  roughness  on  the  fluid  dynamics 
and  heat  transfer  in  the  turbulent  boundary  layer  has  been  in  progress  at 
Stanford  for  the  past  several  years.  Results  of  this  investigation  for  zero 
pressure  gradient  flows  have  been  reported  previously  [1,2].  The  present  study 
considers  the  effects  of  acceleration  on  a turbulent  boundary  layer  in  the 
fully  rough  state.  This  subject  was  investigated  not  only  because  of  its 
importance  in  the  flow  in  nozzles  and  over  turbine  blades  and  reentry 
vehicles,  but  also  to  provide  more  information  on  the  nature  of  turbulence 
by  observing  the  response  of  the  turbulence  field  to  the  imposed  perturba- 
tions (roughness  and  acceleration). 
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1 . 1 Background 

Discussions  of  the  literature  on  turbulent  flows  over  rough  surfaces 
have  been  presented  previously  by  Healzer  [1]  and  Pimenta  [2]  and  will  not 
be  repeated  here.  In  this  section  a brief  introduction  on  the  effects  of 
roughness  on  a turbulent  boundary  layer  will  be  made,  followed  by  brief 
reviews  of  the  zero  pressure  gradient  results  reported  previously  for  the 
present  rough  surface,  the  results  of  accelerated  smooth  wall  turbulent 
boundary  layer  studies,  and  the  few  previous  investigations  of  accelerated 
turbulent  flow  over  rough  surfaces. 

The  influence  of  surface  roughness  on  turbulent  flows  is  usually 
divided  into  three  regimes,  which  are  characterized  by  the  magnitude  of 
the  "roughness  Reynolds  number,"  Re^,  where 


Re 


k 


k U. 


(1.1) 


The  equivalent  sand  grain  roughness  parameter,  k?,  is  a commonly  used,  single- 
length-scale descriptor  of  rough  surfaces  determined  by  comparison  with 
Nikuradse's  [3]  classic  rough  pipe  flow  experiments.  For  Re^.  < 5,  the 
roughness  elements  are  contained  entirely  within  the  viscous  sublayer  and 
the  flow  is  termed  "smooth."  For  5 < Re^  < 55-70  some  of  the  elements 
protrude  through  the  sublayer,  and  the  flow  is  called  "transitionally  rough". 
For  Re^  > 55-70  the  viscous  sublayer  is  effectively  destroyed,  and  the  flow 
is  termed  "fully  rough". 

In  general,  skin  friction  coefficients  and  Stanton  numbers  are  greater 
in  a turbulent  boundary  layer  influenced  by  roughness  than  in  a smooth  wall 
layer  at  the  same  Reynolds  number.  This  causes  larger  temperature  and 
velocity  defects  through  the  layer  and  hence  thicker  boundary  layers,  since 
more  freestream  fluid  is  entrained. 
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Experimental  results  for  zero  pressure  gradient  turbulent  boundary 
layers  on  the  present  rough  surface  were  reported  by  Healzer  [1]  and 
Pimenta  [2].  Healzer  constructed  the  present  experimental  apparatus  and 
reported  C^/2  and  St  data  both  with  and  without  blowing  for  several  veloci- 
ties which  included  the  transitionally  rough  and  fully  rough  flow  regimes. 
He  confirmed  that,  for  fully  rough  flow  over  the  present  surface,  both  C^/2 
and  St  were  independent  of  Reynolds  number,  i.e. 


1 


(1.2) 
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Pimenta  [2]  reported  results  of  an  extensive  investigation  of  the  fluid 
dynamics  and  heat  transfer  in  both  transitionally  rough  and  fully  rough  zero 
pressure  gradient  layers  both  with  and  without  blowing.  His  observations  on 
the  fully  rough  state  included: 

(1)  The  effect  of  roughness  on  the  turbulence  field  structure  extends 
over  most  of  the  layer. 

(2)  Blowing  makes  the  layer  behave  as  if  the  surface  has  physically 
larger  roughness  elements. 

(3)  For  very  large  enthalpy  thicknesses,  the  Stanton  number  appears 
to  converge  to  an  asymptotic  value. 

(4)  Reynolds  shear  stress  correlation  coefficients  are  unchanged  from 
the  values  reported  for  smooth  wall  flows. 

The  response  of  smooth  wall  boundary  layers  to  acceleration  is  discussed 

in  the  summary  report  by  Kays  and  Moffat  [4],  Briefly,  smooth  wall  layer 

dU 

accelerations  are  characterized  by  the  acceleration  parameter,  K = 
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Above  a certain  value  of  K,  the  turbulent  layer  develops  toward  a state 

resembling  laminar  flow.  For  a given  Reynolds  number,  Stanton  number 

2 2 

decreases  with  increasing  K,  and  the  profiles  of  u'  /U  are  lowered  with 
acceleration  [5].  In  a constant  K flow,  the  smooth  wall  turbulent  boundary 
layer  reaches  an  asymptotic  state  where  mean  profiles  are  similar,  Re^  and 
C^/2  are  constant,  and  boundary  layer  thickness  decreases. 

Previously  published  studies  of  the  combined  effects  of  acceleration 
and  roughness  on  the  turbulent  boundary  layer  have  reported  only  values  of 
wall  heat  flux.  Reshocko,  et  al.  [6],  and  Banerian  and  McKillop  [7]  invest- 
igated nozzle  wall  flows,  while  Chen  [8]  cited  experimental  results  for  flow 
over  hemispheres.  No  boundary  layer  information  was  obtained  in  any  of 
these  studies. 

1.2  Objectives 

This  investigation  was  undertaken  to  determine  the  effects  of  accelera- 
tion on  the  fluid  dynamics  and  heat  transfer  in  the  fully  rough  turbulent 
boundary  layer.  Specific  objectives  were: 

(1)  To  define  and  experimentally  verify  the  conditions  required  for 
equilibrium  in  the  fully  rough  turbulent  boundary  layer  with 
pressure  gradient  and  transpiration. 

(2)  To  obtain  a comprehensive  fluid  dynamic  and  thermal  data  set  for 
both  equilibrium  and  nonequilibrium  accelerations  of  the  fully 
rough  turbulent  boundary  layer. 

(3)  To  examine  the  behavior  of  the  mean  and  turbulence  fields  in  the 
accelerated  fully  rough  turbulent  boundary  layer  . 

(4)  To  investigate  the  effect  of  blowing  on  the  equilibrium 
accelerated  layer. 
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1.3  The  Experiment 

A brief  description  of  the  experimental  apparatus  and  measurement 
techniques  will  be  given  in  this  section.  Additional  information  is  provided 
in  Appendix  I . 

The  Stanford  Roughness  Rig  (Figure  1-1)  is  a closed-loop  wind  tunnel 
using  air  as  both  the  primary  and  transpiration  fluids.  Air  temperature  is 
controlled  using  water-cooled  heat  exchangers  in  both  the  primary  and 
transpiration  loops.  The  eight-foot  long,  20-inch  wide  test  section  is  four 
inches  high  at  its  entrance.  A flexible  plexiglass  upper  wall  (constructed 
in  five  sections  connected  by  thin  plexiglass  joints)  can  be  adjusted  to 
give  the  desired  variation  in  U . 

The  test  surface  consists  of  24  plates  each  four  inches  in  the  axial 
direction.  The  plates  (Figure  1-2)  are  0.5  inch  thick  and  uniformly  porous. 
They  are  constructed  of  11  layers  of  0.050- inch  diameter  Oxygen-Free  High 
Conductivity  (0FHC)  copper  spheres  packed  in  the  most  dense  array  and 
brazed  together.  This  configuration  produces  a rough  test  surface  which  is 
uniform  and  deterministic. 

Each  plate  has  individual  electrical  power  and  transpiration  air  controls 
and  thermocouples  for  determining  plate  temperature.  Stanton  number  is 
determined  by  subtracting  the  plate  losses  (known  from  energy  balance 
qualification  tests)  from  the  measured  power  input.  Uncertainty  of  the  St 
data  is  within  + 0.0001  Stanton  number  units  (i.e.,  if  St  = 0.00200,  the 
uncertainty  is  within  + 5%) . 

The  Stanton  number  data  reported  here  were  taken  with  a wall-to- 
freestream  temperature  difference  of  approximately  30°F  to  maintain  a 
constant  property  boundary  layer.  Unless  specifically  stated  otherwise, 
all  St  data  presented  are  for  constant  wall  temperature.  The  freestream 


velocity  at  the  test  section  inlet  was  a nominal  88  ft/sec,  and  all  data 
were  taken  with  a 1/2"  wide,  1/32"  high  phenolic  trip  installed  three  inches 
inside  the  nozzle  exit.  The  turbulent  boundary  layer  was  in  a fully  rough 
state  for  all  cases  reported. 

Mean  temperature  profiles  were  measured  with  a 0.003- inch  diameter, 
butt-welded,  Chromel-constantan  thermocouple  mounted  in  a traversing  probe 
holder.  The  design  was  similar  to  that  of  Blackwell  [9]. 

All  velocity  measurements  were  made  in  an  isothermal  flow  using 
linearized,  constant  temperature  hot-wire  anemometry.  Measurements  of  U 

— 2 — I 

and  u'  were  obtained  using  a horizontal  wire,  while  measurements  of  v'  , 

2 

w'  and  u'v'  were  made  with  a rotatable,  45°  slant  wire. 

The  physical  size  of  the  Roughness  Rig  and  the  porosity  of  the  plates 
imposed  limitations  on  the  strengths  of  the  accelerations  which  could  be 
investigated.  The  height  of  the  tunnel  (four  inches  at  the  nozzle  exit) 
limited  both  the  length  and  severity  of  the  acceleration  region  since 
interference  of  the  top  wall  boundary  layer  with  that  on  the  test  surface 
was  carefully  avoided.  Also,  since  the  plates  were  porous,  the  pressure 
gradient  in  the  axial  direction  induced  flow  through  the  plates  even  with 
the  transpiration  supply  valves  closed.  An  analysis  and  discussion  of  this 
effect  is  presented  in  Appendix  III.  No  effects  of  the  induced  transpiration 
were  apparent  in  the  data.  It  was  concluded  that  the  quantitative  effect  of 
the  induced  transpiration  was  negligible,  certainly  for  the  mildest  and  also 
the  blown  acceleration  runs,  and  that  the  qualitative  trends  in  all  the  data 
(and  the  conclusions  drawn  from  them)  were  unaffected. 
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1.4  General  Organization 


The  general  organization  of  the  results  presented  in  the  following 
chapters  is  described  below.  In  Chapter  2 the  concept  of  "equilibrium"  in 
turbulent  boundary  layers  is  discussed,  and  the  requirements  for  establishing 
equilibrium  in  the  fully  rough  turbulent  boundary  layer  with  pressure 
gradient  and  transpiration  are  developed.  The  experimental  data  are 
presented  in  Chapter  3,  and  characteristics  and  trends  are  discussed.  An 
integrated  discussion  of  the  effects  of  acceleration  on  the  fully  rough 
turbulent  boundary  layer  is  given  in  Chapter  4,  and  Chapter  5 contains  the 
conclusions  of  the  study. 

Additional  information  and  tabular  data  listings  are  contained  in 
Appendices  I - IV. 
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Schematic  of  the  Experimental  Apparatus 


m 4#  €*>  *»■  **•##•  4*4* 

!»••  «*  *>■  ♦*  4k>#>-4 


CHAPTER  2 

THE  EQUILIBRIUM  FULLY  ROUGH  TURBULENT  BOUNDARY 
IAYER  WITH  PRESSURE  GRADIENT  AND  TRANSPIRATION 


The  definition  of  conditions  for  which  a turbulent  boundary  layer 
becomes  similar  in  the  flow  direction  in  some  non-dimensional  sense  has  been 
a subject  of  interest  for  a number  of  years.  Such  similar  behavior  is 
usually  termed  an  "equilibrium"  flow  in  the  literature.  The  term 
"equilibrium"  flow  is  sometimes  used  in  the  sense  of  similarity  of  properly 
normalized  mean  velocity  profiles;  however,  truly  equilibrium  turbulent  flow 
exhibits  similarity  not  only  in  mean  profiles,  but  also  in  the  turbulence 
quantities  themselves. 

The  analytical  and  experimental  work  in  equilibrium  layers  has  been 
prompted  in  part  by  the  desire  to  follow  the  systematic  method  of  fixing  as 
many  variables  as  possible  in  a given  problem.  This  allows  one  to  obtain  a 
better  understanding  of  the  sensitivity  of  the  problem  to  the  uncontrolled 
variables.  In  the  specific  case  addressed  here,  that  of  the  fully  rough 
turbulent  boundary  layer,  the  approach  described  above  is  the  logical  one  to 
follow.  If  the  behavior  of  the  fully  rough  turbulent  boundary  layer 
subjected  to  pressure  gradient  and  transpiration  can  be  examined  in  the 
equilibrium  case,  one  can  then  proceed  to  an  examination  of  the  more 
realistic  non-equilibrium  cases  with  more  confidence  and  understanding. 

In  this  chapter  the  previous  work  in  describing  and  establishing  the 

conditions  for  which  equilibrium  exists  in  turbulent  boundary  layers  is 

10 


discussed.  Following  this,  an  analysis  of  the  fully  rough  layer  is 
presented.  The  analysis  yields  a definition  of  the  proper  acceleration 
parameter  for  fully  rough  flows  and  a description  of  the  conditions  required 
to  establish  equilibrium  in  the  fully  rough  turbulent  boundary  layer  with 
pressure  gradient  and  transpiration. 


2.1  Previous  Studies  of  Equilibrium  Turbulent  Boundary  Layers 


In  1950  Rotta  [10]  presented  an  examination  of  the  conditions  which 
would  yield  a smooth  wall  turbulent  boundary  layer  in  which  the  velocity 
profile  is  distorted  only  affinely  in  the  flow  direction.  He  termed  such 
flows  similar  and  showed  that,  neglecting  the  viscous  wall  region,  the 
equations  describing  the  flow  become  ordinary  differential  equations  if 
C^/2  = constant  and  = axm,  where  a and  m are  constants.  For  a layer 
where  the  friction  coefficient  is  almost  independent  of  x,  similar  solutions 
exist  which  depend  only  on  m and  C^/2 , and  the  boundary  layer  thickness 
increases  linearly  with  x. 

In  1954  Clauser  [11]  presented  experimental  verification  of  the  existence 
of  similar  turbulent  boundary  layer  flows  on  smooth  walls  for  two  different 
adverse  pressure  gradients.  He  termed  such  behavior  "equilibrium”  and 


defined  it  as  the  case  where  both 


R _ 6'  dP 


were  independent  of  x.  In  a later  paper,  Clauser  [12]  showed  that  the 
correct  choice  of  the  length  scale  6'  was  the  displacement  thickness,  so  that 
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It  should  be  noted  that  Clauser's  definition  of  the  shape  parameter  G is 
identical  with  that  presented  earlier  by  Rotta  [10,13].  Additional  theoretical 
treatments  of  equilibrium  turbulent  boundary  layers  were  presented  by 
Townsend  [14]  and  Coles  [15]. 

The  most  definitive  work  on  equilibrium  turbulent  boundary  layer  flow 
was  presented  by  Rotta  [16]  in  1962.  He  showed  that  the  conditions  required 
for  exact  equilibrium  behavior  (reduction  of  the  equations  of  motion  to  an 
ordinary  differential  equation)  are: 


and 


C^/2  = constant 


dfij 
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= constant 


1 dP 

B = — t—  = constant 
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(2.3) 


Two  flows  obeying  these  constraints  exactly  were  shown  to  be  flow  over  a 

smooth  wall  with  U ~ — where  x > x and  flow  over  a uniformly  rough 

00  X _ “ X o 
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lj- I . Other  variations  of  were  shown  to  either 
require  a given  roughness  variation  with  x or  not  to  satisfy  exactly  the 
conditions  required  above. 

There  are  indications  based  on  experimental  rough  wall  studies  that 
exact  equilibrium  cases  exist  for  conditions  not  corresponding  to  the 
velocity  and  roughness  criteria  above.  Perry,  et  al.  [17],  found  that  a 
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wall  with  ~ exp 


zero  pressure  gradient  turbulent  boundary  layer  developing  over  a two- 
dimensional  cavity  type  roughness  of  constant  height  conformed  to  Rotta's 
conditions  for  precise  self-preserving  flow.  Pimenta  [2]  also  found 
indications  in  his  work  on  zero  pressure  gradient  flow  over  the  rough  surface 
used  in  this  study  that  the  boundary  layer  was  approaching  such  an  equilibrium 
state. 

2.2  Conditions  for  Pquilibrium  in  the  Fully  Rough  Turbulent  Boundary  Layer 
with  Pressure  Gradient  and  Transpiration 

In  order  to  determine  the  conditions  for  which  equilibrium  will  be 
obtained  in  the  fully  rough  turbulent  boundary  layer  with  pressure  gradient 
and  transpiration,  consider  the  two-dimens 'onal  momentum  integral  equation 

Cf  dfi..  ! dU 

T+F  = ^r+V2  + H)[r-ar  <2-4) 

oo 

where  the  variation  of  with  x has  been  neglected,  as  have  the  normal 
Reynolds  stresses.  For  the  zero  pressure  gradient  fully  rough  state,  it 
has  been  shown  [1,2]  that  the  skin  friction  is  independent  of  Reynolds 
number  and  can  be  functionally  represented  as 

cf/2  = f {~r>  F)  (2-5) 

where  r is  a length  scale  characteristic  of  the  roughness  elements. 

For  the  present  deterministic  rough  wall  where  height  and  distribution 
are  describable  by  a single  length  scale,  r is  taken  as  the  radius  of  the 
spheres  comprising  the  surface.  In  the  most  general  case,  of  course,  one 
length  scale  describing  height  and  one  describing  distribution  in  addition 
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to  a parameter  describing  roughness  element  form  are  necessary  for  the 
description  of  a rough  surface.  Most  investigators  in  the  past  have  used 
the  "equivalent  sand  grain  roughness"  scale,  ks>  determined  by  comparison 
with  Nikuradse's  [3]  classic  pipe  flow  experiments,  in  order  to  obtain  a 
single  length  scale  description  of  roughness. 

One  condition  necessary  for  equilibrium  is  that  C^/2  be  constant. 
Additionally,  consider  only  the  case  for  constant  F and  assume  that  the 
functional  form  of  Equation  (2.5)  will  remain  valid  for  flows  with  pressure 
gradient.  Under  these  conditions,  62  is  constant  and  Equation  (2.4)  becomes 

6 ^ dUm 

C^/2  + F = (2  + H)  jj-  = constant  (2.6) 


Defining  a pressure  gradient  parameter  for  fully  rough  flow  as 


K 


. dU 
_ _L_  °° 

r " U_  dx 


(2.7) 


where  L is  a length  scale  yet  to  be  specified,  Equation  (2.6)  can  be 
written  as 

(Cf/2  + F) 

Kr  = TTVh7(62/L)  = constant  (2 -8) 

for  an  equilibrium  condition. 

The  choice  of  the  proper  length  scale  L to  use  in  (2.7)  is  not 
immediately  obvious.  One  might  use  an  integral  scale  of  the  flow  (6,  6j,  62) 
or  a roughness  scale  (r,  k^) . The  roughness  element  radius,  r,  will  be  used 
in  this  development.  A discussion  of  the  arguments  for  this  choice  will  be 
deferred  to  a later  section.  Thus, 
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K = _L  00 

r ~ U~  "3x” 


(2.9) 


For  a fully  rough  flow  with  constant  F and  Kr,  the  layer  could  be 
expected  to  exhibit  an  equilibrium  state  for  which  C^/2,  6.,,  H,  and  6 are 
all  independent  of  x.  This  expectation  has  been  experimentally  verified  in 
the  present  investigation  for  positive  Kr  and  F.  For  Kr  < 0 (adverse 
pressure  gradients)  Equation  (2.8)  indicates  equilibrium  flow  is  possible 
only  for  F < 0 (suction).  Fully  rough  flows  with  Kr  constant  are  equilibrium 
flows  in  the  strictest  sense  since  all  of  the  conditions  of  (2.3)  are 
satisfied. 

The  freestream  velocity  variation  required  for  an  equilibrium  flow  is 
found  by  integration  of  Equation  (2.9)  with  = constant  to  be 


U K (x-x  )/r 
00  r o 
1 = e 


(2.10) 


where  the  subscript  o indicates  the  position  at  which  the  velocity  variation 
begins.  This  agrees  with  Rotta's  [16]  result,  but  from  the  development  above 
it  is  clear  that  fully  rough  flow  is  required  for  the  velocity  variation 
(2.10)  to  give  an  equilibrium  flow.  For  transitionally  rough  flow,  C^/2 
is  a function  not  only  of  62/r  and  F,  but  also  of  U^.  Thus, a constant  Kr 
flow  would  not  be  an  equilibrium  flow  for  a transitionally  rough  turbulent 
boundary  layer. 

For  F and  K constant,  it  can  also  be  shown  that 
r 


/ H WV2  + F\ 
~ ' [ITT!  j \-Z^T- ) 


(2.11) 


IS 


Cf/2  = (II  ♦ 2)  Kr(62/r)  - F 


(2.12) 


The  definition  of  Kr  for  fully  rough  flows  is  analogous  to  that  of  the 
smooth  wall  acceleration  parameter 


^2  dx 


(2.13) 


.An  accelerating  turbulent  flow  on  a smooth  wall  with  K = constant  yields  a 

boundary  laver  with  Re.  constant  that  is  equilibrium  in  the  sense  that  mean 

°2 

velocity  profiles  become  similar  and  G and  6 are  constant.  Such  a flow  is 
not  truly  an  equilibrium  flow  in  the  sense  of  equations  (2.3)  since 
d<5.  i 

-r-d-  ~ — j f constant.  A comparison  of  the  asymptotic  accelerated  states  for 
dx  lJz 

oo 

smooth  wall  and  fully  rough  turbulent  boundary  layers  is  presented  in 
Table  2.1. 


2.3  Choice  of  Length  Scale  in  Acceleration  Parameter  K 


The  choice  of  the  correct  length  scale  to  be  used  in  the  fully  rough 
acceleration  parameter  Kr  is  not  obvious  from  the  development  in  Section  2.2. 
A scale  based  on  roughness  size  (r,  k^)  or  a local  scale  of  the  boundary 
layer  (6,  6^,  62)  could  be  chosen.  The  near  wall  scale  used  in  smooth  wall 
layers,  , should  not  be  considered  because  the  turbulence  field  of  the 

T 

fully  rough  layer  is  independent  of  viscous  effects,  at  least  for  regions 
outside  the  roughness  elements  [2]. 

One  requirement  which  should  be  imposed  is  that  when  Kr  is  constant,  an 
equilibrium  condition  should  result.  This  requirement  leads  to  the  choice 
of  roughness  size  as  the  proposed  scale.  Define: 


(2.9) 


K - JL  00 

\ U_  1ST 


6.  dU 

1 °° 

Kr  ~ oi  -ar 


(2.14) 


6,  dU 

v*  _ 2 ® 

Kr  ~ (TUT 


(2.15) 


for  convenience  in  the  discussion  to  follow.  If  the  fully  rough  flow  is  in 

an  equilibrium  state,  6^  and  62  are  both  constant  and  thus  Kr,  K^'.,  and  K*  are 

all  constant  and  meet  the  requirement  above.  However,  consider  a case  where 

a non -equilibrium  acceleration  is  imposed  on  a surface  of  constant  roughness. 

It  is  possible,  in  principle,  that  an  acceleration  could  be  imposed  such  that 

/ 1 dU~\ 

the  product  of  I g — 1 and  6^  or  62  would  be  constant.  Thus,  in  principle, 
Kj.  or  K*  could  be  maintained  constant  in  a non -equilibrium  fully  rough  flow. 
Therefore,  it  appears  that  a local  scale  of  the  layer  is  not  suitable  for 
use  in  defining  K . 

In  choosing  a roughness  length  scale  for  use  in  K , one  is  assuming 

/ 1 <\ 

that  if  r (or  kg)  is  doubled,  then  I y — must  be  halved  to  achieve  the 


must  be  halved  to  achieve  the 


same  effect  for  both  the  cases  r = r^  and  r = 2r^.  Confirmation  of  this 

behavior  must  await  further  experimental  work.  However,  it  is  obvious  that 

some  wall  scale  effect  must  be  included  in  K , otherwise  the  smooth  wall 

r 

parameter,  K,  would  adequately  describe  rough  wall  accelerations. 

Since  the  equivalent  sand  grain  roughness  of  the  present  rough  surface 
according  to  Schlichting  [18]  is  1.25  r,  the  conversion  of  the  values 
reported  to  values  based  on  ks  is  easily  made  if  desired. 
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Table  2.1 


COMPARISON  OF  ASYMPTOTIC  ACCELERATED  STATES  FOR  SMOOTH 
/AND  FULLY  ROUGH  TURBULENT  BOUNDARY  LAYERS 
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Note:  Subscript  o indicates  point  where  acceleration  begins. 
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CHAPTER  3 
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EXPERIMENTAL  RESULTS 

The  experimental  data  obtained  in  this  study  will  be  presented  in  this 
chapter  and  trends,  similarities,  and  comparisons  will  be  noted  and  discussed 
briefly.  An  integrated  discussion  and  description  of  the  effects  of  accelera- 
tion on  the  fully  rough  turbulent  boundary  layer  will  be  presented  in  Chapter 
4. 

The  experimental  program  covered  five  different  cases: 


(1) 

7* 

n 

o 

F = 0 

(2) 

K = 0. 15x10*3 
r 

F = 0 

(equilibrium) 

(3) 

Kr  = 0.29x10*^ 

F * 0 

(equilibrium) 

(4) 

K = 0.29x10*3 
r 

F = 0.0039 

(equilibrium) 

(5) 

K = 0 . 28xl0*6 

F = 0 

(non - equ i 1 i br ium ) 

Case  1 was  run  as  a baseline  set  and  to  compare  the  present  data  with  those 
of  Pimenta  [2]  for  identical  conditions.  Cases  2,  3,  and  4 are  equilibrium 
acceleration  runs  for  the  fully  rough  turbulent  boundary  layer.  In  Case  S 

V dUoo 

the  smooth  wall  acceleration  parameter  K = — j was  maintained  constant. 

00 

This  represents  a non-equilibrium  run  for  the  fully  rough  layer. 

In  setting  up  each  of  the  equilibrium  runs,  the  value  of  Kr  and  the  x- 
position  at  which  the  acceleration  was  begun  were  matched  with  the  H, 

I 
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and  Cj-/2  data  taken  at  that  position  for  = 0,  using  Equation  (2.8).  Thus, 
the  boundary  layer  entered  the  region  of  acceleration  near  the  equilibrium 
state  for  the  Kr  applied,  and  the  length  of  the  equilibrium  flow  established 
was  maximized. 

— 1 — I 

Measurements  included  Stanton  numbers  and  profiles  of  T,  U,  u'  , v'  , 

— 2 

w'  , u'v'.  These  data  allowed  calculation  of  skin  friction  coefficient 
C^/2,  turbulent  Prandtl  number  Pr^.,  mixing  length  l,  and  profiles  of  Q+  and 
t+.  The  profile  measurements  were  obtained  using  two  hot  wires- -one  horizontal 
and  one  rotatable  45°  slant- -and  a butt-welded  thermocouple  probe  similar  in 
design  to  that  used  by  Blackwell  [9],  Details  of  the  measurements  and 
techniques  used  are  presented  in  Appendix  I. 

In  the  following  sections  the  Kr  = 0 baseline  data  will  be  presented 
first.  The  four  cases  with  acceleration  will  then  be  described  with 
presentation  of  the  data  in  the  following  order: 

• Summary  graphs  for  each  case 

• Integral  quantities  ( St,  C^/2,  6,  A,  etc.) 

• Mean  velocity  and  temperature  profiles 

• Reynolds  stress  tensor  components 

• Turbulent  Prandtl  number  and  related  quantities 

The  final  section  of  this  chapter  will  describe  a Stanton  number  prediction 
technique  and  some  supplementary  Stanton  number  data,  including  cases  with 
steps  in  wall  temperature,  variable  wall  temperature,  and  variable  blowing 
with  acceleration. 


3.1  Zero  Pressure  Gradient  Data 

The  data  for  zero  pressure  gradient  were  obtained  both  to  provide  a 
baseline  set  of  measurements  taken  using  the  same  techniques  used  in 
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acquiring  the  accelerated  data  and  also  to  demonstrate  the  compatibility  of 
the  data  with  the  results  of  Pimenta  [2]  for  the  same  conditions  and  test 
surface. 


Stanton  number  data  are  shown  in  Figure  3-1  for  both  F = 0 and  0.0039. 

The  data  of  Pimenta  foi  an  untripped  layer  are  also  plotted,  and  the 
comparison  between  the  two  sets  is  excellent  for  large  — , being  well 
within  the  data  uncertainty  of  + 0.0001  Stanton  number  units.  The  correla- 
tions proposed  by  Pimenta  for  interpolation  of  his  data  are  also  shown.  These 
correlations  are: 


-0.175 


St  = 0.00317 


(3.1) 


for  F = 0 and  4 < — < 15  and 
r 


lnCl+Bb) 


1.175 


d+Bh) 


0.175 


for  0 < F < 0.0040,  where: 


St  is  the  Stanton  number  for  F = 0 and  the  same  A, 
o 2 


= F/St  is  the  blowing  parameter. 

Figure  3-2  presents  the  skin  friction  coefficients  obtained  for  F = 0 

by  Healzer  [1],  Pimenta  [2],  and  the  present  author.  Healzer  differentiated 

his  momentum  thickness  measurements  to  obtain  C^/2,  while  Andersen's  [19] 

shear  stress  method  for  skin  friction  determination  was  used  in  this  study 

and  also  by  Pimenta.  The  results  of  Pimenta  and  the  present  author  show 

good  agreement,  while  the  data  of  Healzer  deviate  slightly  from  the  others 

at  the  larger  values  of  67/r.  The  correlation  of  Pimenta  for  F = 0 and 
62 

1.0  < — < 10.0 

21 


Cf/2  = 0.00328  (y) 


(3.3) 


-0.175 


is  also  plotted. 

All  of  the  skin  friction  coefficients  in  this  study  were  calculated 

using 


(3.4) 


1: 


The  derivation  of  (3.4)  is  straightforward.  Briefly,  the  momentum 
equation  (incorporating  the  usual  boundary  layer  assumptions  but  allowing 
pro  = p^fx))  and  the  continuity  equation  are  integrated  from  the  surface  to  a 
position  y^  in  the  boundary  layer.  Then,  measurement  of  successive  velocity 
profiles  in  the  x-direction  and  u'v'  at  y = y^  for  each  x-position  allows 
calculation  of  C^/2  versus  x,  using  (3.4).  The  position  y^  was  always  0.130" 
in  this  study,  since  the  rotatable  slant  hot  wire  used  to  measure  u'v'  was 
limited  to  y > 0.125". 

A typical  velocity  profile  is  plotted  versus  (y  + Ay)/^  in  Figure  3.3. 
Since  the  normal  coordinate  y is  referred  to  the  tops  of  the  spherical 
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elements  comprising  the  test  surface,  the  "wall  shift"  Ay  gives  the  location 
of  the  "apparent  wall"  for  the  mean  velocity.  This  wall  shift  has  been  a 
topic  of  much  discussion  by  previous  workers  in  roughness  and  will  be 
discussed  in  more  detail  later  in  this  chapter.  It  is  shown  in  the  figure 
that  the  present  data  follow  Schl icht ing ' s (18]  expression  for  fully  rough 
flow 

oL  ■ i? ln  f * 8-5 

T S 

Note  that  the  value  of  k^  used  (0.031")  is  determined  from  Schlichting's 

tabulated  values  and  not  by  back-fitting  Equation  (3.5).  The  smooth  wall 

"law  of  the  wall"  is  also  shown  for  reference. 

Measurements  of  the  three  components  of  the  turbulent  kinetic  energy 

normalized  by  IT  are  shown  in  Figure  3-4  plotted  versus  y/< 5.  The  present 

measurements  agree  with  those  of  Pimenta  within  the  data  uncertainty. 

2 

Comparison  of  the  fully  rough  data  with  the  u'  data  of  Klebanoff  [20]  for 
a smooth  wall  shows  several  important  characteristics  of  fully  rough  flow 

(which  were  noted  by  Pimenta  [2]).  First,  for  fully  rough  flow  the  peak  in 

2 

u'  is  moved  out  from  the  wall,  lowered,  and  spread  over  a greater  portion 
of  the  layer  than  is  the  case  for  smooth  wall  flows.  Second,  the  effect  of 
the  roughness  is  felt  across  practically  the  entire  layer  in  the  form  of 
increased  turbulence  energy.  Blake  [21]  also  observed  this  behavior  in  his 
fully  rough  flow  data.  Thus,  the  assumption  made  by  some  authors  [8,17]  that 
the  effect  of  roughness  is  confined  to  the  near  wall  region  is  not  valid  for 
the  turbulent  kinetic  energy  components.  Pimenta  showed  that  this  effect 
was  not  due  to  the  greater  freestream  turbulence  (=  0.41)  in  the  present 
tunnel  as  compared  with  that  of  Klebanoff  (~  0.02%).  He  also  showed  that 


i 


a 


I 

I 
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the  use  of  as  a normalizing  velocity  did  not  collapse  the  smooth  and 
rough  wall  results,  as  was  suggested  by  Hinze  [22]  based  on  the  measurements 
of  Corrsin  and  Kistler  [23]  over  2-D  roughness  elements. 

3.2  Data  with  Acceleration 

Summary  graphs  for  the  four  cases  of  accelerated  flow  investigated  are 
shown  in  Figures  3-5  through  3-8.  The  purpose  of  these  plots  is  to  show 
the  variation  of  Kr  and  the  integral  quantities  C^/2,  62,  and  H which  are 
indicators  of  equilibrium  flow  according  to  the  discussion  in  Chapter  2.  The 
Stanton  number  variation  is  also  plotted  to  illustrate  the  integral  behavior 
of  the  thermal  field.  In  each  figure,  he  data  are  plotted  versus  distance 
along  the  test  section,  x.  In  the  discussion  which  follows,  F = 0 unless 
specifically  stated  otherwise. 

Data  from  the  Kr  = 0.15  x 10  ^ run  are  presented  in  Figure  3-5.  This 
run  had  the  longest  region  of  Kr  constant  (x  = 44”  to  88”),  with  the  velocity 
increasing  from  approximately  88  to  115  ft/sec.  As  seen  in  the  figure,  62, 

H and  C^/2  all  become  constant  in  the  region  of  constant,  indicating  that 
equilibrium  flow  was  established.  Stanton  numbers  in  the  acceleration  region 
are  about  10%  larger  than  for  the  Kr  = 0 case  and  appear  to  be  approximately 
constant  within  the  data  uncertainty.  The  behavior  is  different  from  that 
observed  for  accelerated  smooth  wall  layers,  where  Stanton  number  is 
unaffected  for  small  K,  then  decreases  in  comparison  to  the  unaccelerated 
case  at  the  same  Reynolds  number  or  same  x-position  as  K increases  14,24,25,26]. 

Data  for  = .29  x 10  ^ are  shown  in  Figures  3-6  and  3-7  for  F = 0 
and  F = 0.0039,  respectively.  In  both  these  cases,  Kr  is  constant  from 
x = 24"  to  52”,  increases  from  88  to  129  ft/sec,  and  ^ C^/2,  and  H all 
reach  constant  values  in  the  acceleration  region.  Stanton  number  shows  the 
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same  behavior  as  seen  in  Figure  3-5  in  the  region  of  acceleration,  then 
decreases  immediately  to  the  K}.  = 0 baseline  data  when  the  acceleration  is 
removed . 

The  summary  data  for  the  K = .28  x 10  6 , nonequilibrium  case  are 
presented  in  Figure  3-8.  The  smooth  wall  acceleration  parameter  K is 
constant  from  x = 24"  to  52",  U increases  from  88  to  150  ft/sec,  and  K 

'co  ' Y 
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varies  from  .25  - .50  x 10  in  this  region.  The  shape  parameter  K 
decreases  along  the  entire  test  section,  while  62  increases  as  the  layer  is 
entering  the  region  of  acceleration,  then  levels  off  and  finally  decreases. 
This  &2  behavior  is  similar  to  that  observed  in  the  asymptotic  accelerated 
smooth  wall  layer  [27,28].  Skin  friction  coefficient  shows  very  little 
variation,  and  appears  to  remain  about  constant.  This  is  not  surprising 
considering  the  small  variation  of  62  in  the  acceleration  region.  Stanton 
number  shows  the  same  increase  over  Kr  = 0 values  observed  in  the  equilibrium 
cases  and  recovers  immediately  to  unaccelerated  baseline  values  when  the 
acceleration  is  removed. 

It  was  noted  above  that  Stanton  number  appears  to  be  approximately 
constant  within  the  data  uncertainty  in  regions  of  Kr  constant.  However, 
it  is  impossible  to  reach  a firm  conclusion  in  this  regard  due  to  the 
inherent  uncertainty  in  the  data  and  the  relatively  short  regions  of 
acceleration.  An  argument  that  Stanton  number  for  Kr  > 0 varies  as  A2  ^’^5 
(as  in  the  K = 0 case)  or  some  similar  weak  function  of  A?  could  also  be 
supported  by  the  present  data. 

If,  for  the  sake  of  argument,  one  assumes  that  St  is  constant  in  a 
region  of  constant  Kr,  then  the  behavior  of  A2  in  such  a region  can  be 
determined.  The  2-D  energy  integral  equation  for  F = 0,  constant  properties, 
and  constant  wall  to  freestream  temperature  difference  may  be  written  as  [29] 


(3.6) 


If  St  and  Kr  are  constant  from  a position  xc  to  x,  then  (3.6)  can  be 
integrated  to  yield 
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Thus  the  enthalpy  thickness  will  approach  an  asymptotic  value  (Stcr/Kr)  if  a 

flow  is  established  such  that  St  and  K are  constant. 

r 
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Equation  (3.7)  was  evaluated  for  the  Kr  = .15  x 10  case  assuming 
xc  = 58"  and  Stc  = 0.00242.  The  results  are  shown  in  Figure  3-9  and 
compared  with  the  enthalpy  thicknesses  computed  from  the  mean  velocity  and 
temperature  profiles.  The  agreement  between  the  measured  values  and 
calculated  values  assuming  St  = constant  is  excellent,  thus  supporting  the 
observation  that  Stanton  number,  if  not  a constant,  is  at  most  a weakly 
varying  function  of  A£  in  an  equilibrium  accelerated  fully  rough  turbulent 
boundary  layer  with  constant  wall  temperature. 

_3 

As  shown  in  Figure  3-9,  for  the  = .15  x 10  run  the  approach  of  b>2 
to  the  indicated  asymptotic  value  is  very  slow.  In  fact,  A2  would  reach  99% 
of  the  asymptotic  value  only  after  (x-xc)  « 60  feet.  Since  62  is  constant 
in  the  equilibrium  case,  the  ratio  ^2^2  wou^  therefore  increase  for  an 
appreciable  distance. 


3.2.1 


The  Stanton  number  data  in  the  accelerated  region  for  the  three 


equilibrium  runs  are  shown  versus  in  Figure  3-10  compared  with 

unaccelerated  data  for  F = 0 and  0.0039.  The  accelerated  data  increase 
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over  the  Kf  = 0 data  by  - 10%  for  F = 0 and  ~ 20%  for  F = 0.0039.  As 
discussed  in  the  previous  section,  the  Stanton  number  data  ir.  the  constant 
Kr  region  can  be  argued  to  be  either  approaching  a constant  value  or  at  the 
very  least  to  be  a weaker  function  of  enthalpy  thickness  than  in  the 
unaccelerated  case. 

Skin  friction  coefficients  in  the  acceleration  region  are  plotted 
versus  Fig1-11"6  3-11  and  compared  with  the  unaccelerated  cases.  The 

Kr  = 0 data  for  F = 0.0039  are  from  Pimenta's  study.  For  the  three  equili- 
brium cases,  both  C^/2  and  are  constant  in  the  acceleration  region  and 
thus  only  a single  data  point  for  each  case  appears  in  these  coordinates. 

In  the  unblown  cases,  it  appears  that  acceleration  causes  a slight  increase 
(~  5%)  in  C^/2  over  the  baseline  data.  It  should  be  noted  that  this  is 
within  the  uncertainty  (~  10%)  of  the  C^/2  data,  however.  In  the  blown  case 
(F  = 0.0039),  the  acceleration  data  point  lies  approximately  30%  above 
Pimenta's  Kr  = 0 data.  In  smooth  wall  flows  acceleration  also  leads  to  an 
increase  in  C^/2  compared  with  zero  pressure  gradient  values  for  the  same 
Reynolds  number  [27,28]. 

Figure  3-12  shows  temperature  and  velocity  boundary  layer  thicknesses 
(A  and  6,  respectively)  for  the  five  cases  investigated.  No  temperature 
profiles  were  taken  in  the  present  study  for  = 0,  so  A is  not  shown  for 
the  baseline  case.  In  the  three  equilibrium  accelerated  cases,  the  rate  of 
growth  of  both  6 and  A decreases  in  the  acceleration  region.  From  the 
equilibrium  conditions  developed  earlier  for  the  fluid  dynamics,  one  would 
expect  6 to  eventually  assume  a constant  value  for  Kr  = constant,  and  it 
does  appear  from  Figure  3-12  that  6 is  approaching  an  asymptote  in  the  region 
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of  constant.  After  the  acceleration  is  removed  in  the  two  Kr  = 0.29  x 10 

runs,  the  boundary  layer  thicknesses  A and  5 resume  a rate  of  increase  with  x 

similar  to  that  observed  for  the  Kr  = 0 case.  The  data  for  the  nonequilibrium 
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K = .28  x 10  ^ run  exhibit  behavior  similar  to  that  observed  in  smooth  wall 
accelerated  flows  [24,25,27,28].  Both  6 and  A begin  to  decrease  near  the 
end  of  the  acceleration  region.  In  all  four  acceleration  cases,  the 
temperature  boundary  layer  thickness  A is  greater  than  6 for  all  x,  but  the 
two  thicknesses  show  the  same  trends  in  the  acceleration  and  recovery  regions. 

A comparison  of  enthalpy  thicknesses  obtained  from  integration  of 
temperature  and  velocity  profiles  and  from  integration  of  the  constant  property 
energy  integral  equation  in  the  form 


St  + F * 


U (T  -T  J dx 
°°v  w °°,o' 


A-  U (T  -T  ) 
2 oo  v w °o,0 


(3.8) 


using  measured  Stanton  numbers  is  shown  in  Figure  3-13  for  the  four  acceler- 
ation cases.  Reasonable  agreement  is  found  between  the  two  methods,  with 
the  maximum  discrepancy  being  about  10%.  The  behavior  of  similar  to 

that  observed  for  A in  Figure  3-12. 

Figure  3-14  presents  the  variation  of  roughness  Reynolds  number  with  x 
for  all  five  cases,  where 


Rek 


kS  UT 


(3.9) 


and  kg  was  taken  as  0.031",  as  noted  previously,  for  the  present  surface. 

The  roughness  Reynolds  number  increases  with  acceleration  since  U (=/Cf/2  UM) 
k. 

increases  and  remains  constant. 

These  results  have  important  implications.  The  utility  of  the  roughness 
Reynolds  number  lies  in  its  magnitude  relative  to  the  viscous  sublayer 
thickness.  Following  the  traditional  argument,  for  Re^  < 5,  the  roughness 
elements  do  not  penetrate  the  sublayer  and  the  flow  retains  its  smooth  wall 
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characteristics.  For  5 < Re^  < 55  to  70  (depending  on  the  data  and/or 
author)  the  flow  is  "transitionally"  rough,  and  for  Re^  > 55  to  70  the  flow 
is  fully  rough.  These  ranges  are  all  for  F = 0.  Since  Re^,  increases  in  the 
acceleration  region,  the  roughness  elements  protrude  further  out  into  the 
layer  (in  a nondimens ional  sense)  in  this  region.  There  is  no  viscous 
sublayer  present  in  the  fully  rough  layer,  so  the  increase  in  Re^  with 
acceleration  can  be  viewed  as  making  it  more  difficult  for  a viscous  sublayer 
to  form. 

This  observation  is  important  when  one  considers  the  behavior  of 
accelerated  smooth  wall  flow.  Kays  and  Moffat  [4]  note  that  experimental 
evidence  indicates  acceleration  of  a smooth  wall  turbulent  boundary  layer 
causes  an  increase  in  the  viscous  sublayer  thickness.  Also,  it  is  well  known 
from  the  results  of  many  investigations  that  acceleration  of  a smooth  wall 
turbulent  layer  causes  the  layer  to  develop  toward  a state  resembling  laminar 
flow.  Consideration  of  these  smooth  wall  accelerated  flow  characteristics 
might  lead  one  to  expect  a fully  rough  turbulent  boundary  layer  subjected  to 
a favorable  pressure  gradient  to  develop  first  transitionally  rough,  then 
finally  smooth  wall  characteristics.  The  present  results  indicate  that  this 
is  not  the  case.  To  the  contrary,  acceleration  makes  a fully  rough  flow 
appear  "rougher"  in  the  sense  that  the  roughness  elements  protrude  further, 
nondimensionally,  into  the  turbulent  layer. 


3.2.2  Mean  Velocity  and  Temperature  Profiles 


Figures  3-15  through  3-18  present  mean  velocity  profiles  for  the 


four  acceleration  runs  plotted  as  (J/U^  versus  y/6 -, . In  these  and  subsequent 
figures,  x denotes  the  x position  at  which  the  relevant  acceleration  para- 
meter (Kr  or  K)  becomes  constant. 
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In  the  graphs  for  the  three  equilibrium  runs  (3-15  through  3-17)  the 
profiles  are  similar,  as  expected,  after  the  layer  is  a sufficient  distance 
into  the  acceleration  region.  The  similarity  extends  down  to  the  first 
point  from  the  surface  (y  = 0.006").  The  nonequilibrium  data  (Figure  3-18) 
do  not  exhibit  similarity. 

Figures  3-19  and  3-20  present  profiles  from  the  K = 0.15  x 10  ^ and 
unblown  Kr  = 0.29  x 10  ^ cases,  respectively,  plotted  as  U/V  vs  (y  + A y)/^- 
The  smooth  wall  "law  of  the  wall"  and  Fchlichting’s  [18]  expression  for 
fully  rough  flow  (Equation  3.5)  are  also  shown  for  comparison.  As  noted 
previously  in  this  chapter,  the  wall  shift  Ay  locates  the  apparent  or  virtual 
location  of  the  surface  below  the  tops  of  the  roughness  elements. 

The  wall  shift  was  determined  by  the  same  technique  used  by  Pimenta  [2], 
i.e.,  the  method  suggested  by  Monin  and  Yaglom  [30].  Briefly,  if  it  is 
assumed  that  a logarithmic  law  of  the  wall  region  exists  in  the  velocity 
profile,  it  can  be  shown  that 


_U_ 

U 


(3.10) 


where 


k = Karman  constant  (*  0.41) 

z = constant 
o 

Ay  = constant 

The  proper  wall  shifts  were  determined  using  a form  of  (3.10)  - Ay  was  varied 
until  a value  was  determined  for  which  z was  constant. 
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It  was  found  that  Ay  = 0.006"  for  all  the  profiles  in  the  present 
unblown  data.  Since  this  is  the  same  value  found  by  Pimenta  for  his  zero 
pressure  gradient  data,  it  can  be  concluded  that,  for  the  Kr  range  of  this 
study,  Ay  is  unaffected  by  favorable  pressure  gradients  and  does  not  vary 
with  x.  This  result  is  quite  different  than  that  reported  by  Perry,  et  al. 
117],  who  investigated  turbulent  boundary  layer  flow  over  2-D  roughness 
elements  for  both  zero  and  adverse  pressure  gradients.  They  found  that  Ay 
varied  with  x,  and  in  fact,  was  actually  larger  than  the  roughness  height 
under  some  adverse  pressure  gradient  conditions. 

Comparison  of  the  present  profiles  with  Schlichting' s expression  shows 

that  the  constant  would  have  to  be  increased  from  8.5  to  approximately  9.1 

to  match  the  accelerated  data.  The  reason  for  this  shift  is  not  known.  The 

decrease  in  the  value  of  AU/U^  between  the  smooth  wall  law  of  the  wall  and 

the  present  data  when  acceleration  is  imposed  should  not,  in  the  author's 

opinion,  be  taken  as  an  indication  the  flow  is  tending  toward  the  transi- 

~~2 

tionally  rough  state.  The  u'  profiles  to  be  presented  later  in  this  chapter 
exhibit  none  of  the  transitionally  rough  characteristics  described  by 
Pimenta  [2]  for  this  surface.  In  addition,  the  increase  of  Re^  in  the 
acceleration  region  indicates  a trend  away  from,  rather  than  toward  the 
transitionally  rough  state  (see  Section  3.2.1). 

A comparison  between  the  blown  and  unblown  velocity  profiles  for 
Kr  = 0.29  x 10  ^ is  presented  in  Figure  3-21.  The  behavior  is  as  expected  - 
the  injection  of  low  momentum  fluid  at  the  wall  lowers  the  mean  velocity 
compared  to  the  unblown  case. 

Temperature  profiles  for  the  three  equilibrium  acceleration  cases 
are  plotted  in  Figures  3-22  through  3-24  in  CTW_T') / (Tw~Too)  vs  y/&2 
coordinates.  Similarity  of  the  profiles  is  observed  in  all  cases  and 
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extends  to  the  closest  data  point  from  the  surface  (y  = 0.013").  Temperature 
profiles  for  the  K = 0.28  x 10  ^ non -equilibrium  run  are  shown  in  Figure  3-25. 
The  apparent  similarity  observed  here  is  surprising  but  can  be  explained  by 
reference  to  the  behavior  of  the  thermal  boundary  layer  thickness  A and 
enthalpy  thickness  for  this  case  shown  in  Figures  3-12  and  3-13,  respec- 
tively. Both  A and  A^  vary  little  in  the  acceleration  region  - they  appear 
to  be  approaching  maxima,  and  a decrease  is  actually  observed  in  the  data 
for  A at  the  end  of  the  acceleration  region.  Thus  the  combination  of 
approximately  constant  A and  t\^  and  approximately  constant  Stanton  number 
(Figure  3-8)  in  the  region  of  acceleration  leads  to  similarity  in  the 
temperature  profiles.  This  similarity  would  probably  not  be  maintained  if 
the  region  of  K constant  were  extended. 

Temperature  profiles  for  the  Kr  = 0.15  x 10  ^ run  are  plotted  in 
Figure  3-26  in  (T^-T) / (Tw-T<x>)  versus  U/U^  coordinates.  Pimenta  [2]  found 
these  coordinates  useful  since  Kr = 0 data  are  linear  when  plotted  in  this 
manner.  The  present  profile  at  x = 34"  (prior  to  the  acceleration  region) 
exhibits  this  linearity.  The  two  profiles  in  the  acceleration  region, 
however,  are  not  linear  and  do  not  exhibit  similarity  in  these  coordinates. 

It  should  be  noted  that  the  accelerated  profiles,  if  extrapolated  to  U/U^  = 0, 
still  show  the  temperature  "jump"  condition  discussed  by  Pimenta,  indicating 
that  the  apparent  wall  position  is  different  for  the  temperature  and  velocity 
fields. 

A comparison  of  the  blown  and  unblown  temperature  profiles  for 
Kr  = .29  x 10  3 is  shown  in  Figure  3-27.  The  results  are  as  expected,  with 
the  injection  of  fluid  at  the  wall  temperature  resulting  in  higher  tempera- 
tures in  the  near  wall  region. 
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3.2.3  Reynolds  Stress  Tensor  Components 


The  behavior  of  the  turbulence  quantities  will  be  examined  in 
this  section.  Due  to  the  physical  limitations  of  the  apparatus,  only  in 
the  Kr  = 0.15  x 10  ^ run  was  an  acceleration  region  of  sufficient  length 
established  to  investigate  the  similarity  of  the  turbulence  quantities  in 
the  flow  direction.  The  data  from  this  run  will  form  the  primary  basis  for 
the  discussion  of  the  effects  of  acceleration  on  the  turbulence  field.  Data 
from  the  other  runs  will  be  presented  as  additional  support  for  the  points 
presented;  however,  a direct  comparison  of  the  exact  magnitude?  between  the 
data  of  different  accelerations  is  not  particularly  meaningful  due  to  the 
variation  in  the  values  of  (x-x  )/5  at  which  the  profiles  of  the  different 

d 

runs  were  taken. 

In  the  following  discussion,  the  similarity  of  turbulence  profiles  in 
the  flow  direction  is  examined  first.  Comparisons  are  then  made  between  the 
acceleration  profiles  and  those  for  Kr  = 0,  and  finally  a comparison  of  the 
Reynolds  shear  stress  correlation  coefficients  for  the  different  runs  is 
presented.  A comparison  of  correlation  coefficients  between  the  different 
acceleration  runs  is  valid  due  to  the  demonstrated  insensitivity  of  the 

coefficients  to  the  boundary  conditions  imposed. 

.3 

Figure  3-28  presents,  for  Kr  = 0.15  x 10  , the  three  components  of  the 

turbulent  kinetic  energy  nondimensionalized  by  IT  and  plotted  versus  y/6 2 
for  (x-x  ) = 22  and  42  inches.  Excellent  similarity  is  observed  in  the  u' 

a 

component,  and  the  agreement  in  the  v'  and  w'  components  is  within  the 
uncertainty  of  the  data  (~  10%).  The  Reynolds  stress  profiles  at  the  two 
positions  are  compared  in  Figure  3-29.  The  difference  between  the  two 
profiles  is  on  the  order  of  the  data  uncertainty.  The  results  shown  in 
Figure  3-28  and  3-29  demonstrate  that  a state  of  similarity  is  being 
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approached  by  the  turbulence  quantities,  with  the  data  indicating  that  u'v' 
possibly  requires  a greater  distance  to  become  truly  similar  than  do  the 
other  quantities. 

2 2 

Profiles  of  u'  /U^  versus  y/<5  are  shown  in  Figure  3-30  for  the  Kr  = 0 

and  .15  x 10  ^ cases.  The  decrease  in  longitudinal  turbulence  intensity  with 

acceleration  is  quite  evident  and  similar  to  the  behavior  observed  with 

accelerated  smooth  wall  flows  [5].  When  the  profiles  in  Figure  3-30  are 

compared  with  the  same  two  profiles  in  Figure  3-31  (where  the  data  are 
2 — I 

normalized  by  U^) , one  observes  that  the  peaks  in  u'  nearly  coincide  when 

2 2 

scaling  is  used  but  are  displaced  in  level  if  U^_  scaling  is  used.  This 

— 1 2 

near -co incidence  of  the  u’  peaks  when  normalized  by  Uro  provides  a convenient 
reference  level  when  canparing  the  profiles,  and  all  fluctuation  data  to 
fol low  are  presented  in  this  form. 

The  three  components  of  the  turbulent  kinetic  energy  for  the  = 0 and 

.3  ? 

0.15  x 10  cases  are  compared  in  Figure  3-31  as  u!  versus  y/6.  As  stated 

above,  the  level  of  the  u'  component  in  these  coordinates  is  changed  very 
little  by  acceleration  for  y/6  < 0.1.  The  v'  and  w'  components  are  sub- 
stantially lower  than  the  Kr  = 0 data  in  the  region  y/6  « 0.1,  while  in  the 
outer  region  (y/6  ;>  0.2)  all  three  components  are  lowered  on  the  order  of  40% 
compared  to  the  Kr  = 0 values.  Thus,  when  compared  with  the  unaccelerated 
data,  acceleration  decreases  the  level  of  turbulent  kinetic  energy  over  the 
entire  layer  and  makes  the  turbulence  structure  much  more  anisotropic  in  the 
inner  region.  Unfortunately,  no  measurements  of  the  v'  and  w'  components  in 
a smooth  wall  accelerated  layer  are  known  to  the  author,  so  no  comparison  of 
rough  and  smooth  wall  behavior  with  acceleration  is  possible. 

Consideration  of  the  turbulent  kinetic  energy  equation  and  the 
equations  for  the  energies  in  the  three  components  allows  some  insight  into 
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the  behavior  observed  in  Figure  3-31.  The  time  averaged  turbulent  kinetic 
energy  equation  for  stationary  flow  and  no  body  forces  can  be  written  as 
[30,31] 


) IE  l).  + i p uTuTuT  + n'u!  - u'.o!  . | 
r-  \ 1 2 J 3 i 1 i J ij  j 
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where 
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The  terms  on  the  left  side  of  (3.11)  are,  respectively,  the  spatial 
transfer  of  E by  the  mean  motion,  by  the  turbulence  fluctuations,  by  the 
"pressure  diffusion",  and  by  the  viscous  shear  stresses  of  the  turbulence 
field.  The  terms  on  the  right  hand  side  are  the  dissipation  of  E by 
molecular  viscosity  and  the  production  of  E by  the  interaction  of  the 
Reynolds  stress  tensor  with  the  mean  velocity  gradients. 

The  equations  [30]  for  the  three  components  of  E contain  terms  similar 
to  those  in  (3.11)  and,  in  addition,  the  terms  p'  p'  and  p' 

— 2 — 2 — I 

appear  on  the  right  hand  side  of  the  u'  , v'  , and  w'  equations,  respec- 
tively. Since  these  three  additional  terms  sum  to  zero  by  continuity,  they 
do  not  appear  in  the  equation  for  the  total  turbulent  kinetic  energy.  Thus 
these  pressure  fluctuation- turbulence  field  interaction  terms  transfer  energy 


among  the  components  of  E,  but  play  no  direct  role  in  the  spatial  transfer 
of  turbulence  energy. 


In  the  flows  of  this  investigation,  all  the  terms  of  the  production 
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u.u. 
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. ....  ~ — are  negligible  except  u'v'  -r—  (see  Section  3.2.4).  Therefore  the 

i j <>y  

entire  turbulent  kinetic  energy  production  goes  into  the  u'  component  of  E, 


T? 


.,.2 


and  the  v'  and  w'  components  receive  energy  only  through  the  pressure 

fluctuation- turbulence  field  interaction  terms  described  above.  Since  the 

effect  of  acceleration  is  to  make  the  fully  rough  layer  much  more  anisotropic 

in  the  inner  region,  acceleration  must  decrease  the  sum  of  the  pressure 

— 1 — 2 

fluctuation  transfer  (source)  and  dissipation  (sink)  terms  in  the  v'  and  w' 
equations.  This  argument  can  be  carried  further  only  if  one  assumes  the 


dissipation  is  affected  only  slightly  by  acceleration- -under  this  assumption, 

3v’ 

9y 


3v ' 3w ' 

it  would  have  to  be  true  that  the  correlations  between  p'  and  are 


decreased  significantly  by  acceleration. 

2 

Profiles  for  the  components  of  q are  presented  in  Figures  3-32  and  3-33 
for  the  Kr  = 0.29  x 10  3,  F = 0,  and  K = 0.28  x 10  6 cases,  respectively. 

The  behavior  observed  in  each  case  is  similar  to  that  already  shown  in 
Figure  3-31  for  the  Kr  = 0.15  x 10  3 data. 

2 

The  effects  of  acceleration  on  the  components  of  q in  a fully  rough 
layer  with  blowing  are  shown  in  Figure  3-34.  The  data  points  are  for  the 
Kr  = 0.29  x 10  3,  F = 0.0039  case  of  the  present  study,  while  the  solid 

lines  for  Kr  = 0,  F = 0.0039  are  from  Pimenta  [2].  Compai  son  of  the  two 

sets  of  data  yields  several  important  points.  First  acceleration  decreases 
all  components  of  q in  the  outer  region  much  as  it  does  in  the  unblown 

layers.  However,  the  behavior  in  the  inner  region  is  quite  different  from 

that  in  the  unblown  cases.  The  degree  of  anisotropy  in  the  accelerated  data 
is  about  the  same  at  y/6  = 0.2  as  in  the  zero  pressure  gradient  data. 


Unfortunately,  the  probe  size  restrictions  prevented  acquisition  of  v'  and 

w'  data  for  y/6  <0.1  and  the  trends  inside  this  region  are  undetermined.  I 

In  the  unblown  cases,  the  general  shapes  of  the  profiles  with  and 
without  acceleration  were  similar.  However,  in  Figure  3-34  one  sees  that  in 
the  blown  layer  acceleration  alters  the  basic  shape  of  the  profiles.  The 
curvature  of  the  u'  profile  in  the  outer  region  is  completely  changed,  and 
the  peaks  in  the  v'*"  and  w'“  profiles  for  Kr  = 0 are  suppressed.  In  fact, 
the  profiles  from  the  blown  accelerated  layer  look  much  more  like  the 

unblown  K = 0 profiles  than  the  profiles  from  the  blown  Kr  = 0 case.  This 

2 a 

can  be  seen  in  Figure  3-35  where  the  u'  component  is  shown. 

2 

The  Reynolds  shear  stress  nondimens ionali zed  by  U is  shown  in  Figure 

3-36  for  the  Kr  = 0 and  0.15  x 10  ^ cases.  The  behavior  of  u' v'  with 

acceleration  is  similar  to  that  calculated  for  the  smooth  wall  accelerated 

layer  [25].  No  measurements  of  this  term  in  the  smooth  wall  accelerated 

layer  have  been  published  to  the  knowledge  of  the  author.  The  observed 
2 

decrease  in  -u'v'/U  with  acceleration  would  lead  one  to  expect  a probable  j 

decrease  in  the  production  of  turbulent  kinetic  energy  with  acceleration  (in 
a nondimensional  sense).  This  point  will  be  expanded  in  the  following 
section. 

Figure  3-37  presents  the  measured  correlation  coefficients  Ruv  and  R 2 

where 

R = - uV*' 
uv 

and 


(3.13) 
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The  measured  values  for  all  four  acceleration  cases  are  in  excellent  agree- 


ment with  those  (Ruy  * 0.45,  ~ 0.15)  reported  for  both  smooth  wall 

layers  [52,33,34]  and  zero  pressure  gradient  rough  wall  layers  [2],  It  thus 
appears  that  the  relationship  between  the  Reynolds  shear  stress  and  the 
diagonal  components  of  the  tensor  is  truly  universal  and  independent  of 
boundary  conditions. 

Since  the  turbulent  shear  and  turbulent  kinetic  energy  are  primarily 
generated  during  periods  of  bursting  [35,36,37],  it  is  logical  to  propose 
that  the  universal  values  of  Ruv  and  R^2  observed  result  from  a universal 
attribute  of  the  bursting  and  decay  process  itself.  Grass  [38],  who  reported 
results  of  a hydrogen  bubble  technique  investigation  of  a turbulent  water 
channel  flow  over  a pebble-type  rough  surface,  observed  that  the  bursting 
process  appeared  more  vigorous  in  the  fully  rough  than  the  smooth  wall  case. 
The  inrushing  fluid  interacted  with  the  fluid  among  the  roughness  elements 
(which  is  more  energetic  than  that  in  the  viscous  sublayer  on  a smooth  wall) , 
and  in  the  ejection  phase  of  the  process  the  fluid  moved  almost  vertically 
upward.  These  results  are  consistent  with  Pimenta's  results  of  higher 
turbulence  energy  throughout  the  layer  in  the  fully  rough  state.  Thus,  the 
levels  of  shear  stress  and  energy  are  influenced  by  the  vigor  of  the  bursting 
process  and,  by  extension,  the  boundary  conditions.  However,  it  appears  from 
all  the  data  available  that  in  any  flow  where  the  level  of  turbulence  is 
generated  and  maintained  by  the  bursting  process,  the  relationship  between 
the  components  of  the  Reynolds  stress  tensor  is  fixed  by  some  basic  attribute 
of  the  bursting  and  decay  mechanisms  and  is  independent  of  boundary  conditions. 
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3.2.4  Turbulent  Prandtl  Number  and  Related  Quantities 


I 

I 

I 


The  results  discussed  in  this  section  were  obtained  from  calcu- 
lations using  St,  C^/2,  l),  and  T data  and  the  energy,  momentum,  and  continuity 
equations  integrated  to  a position  in  the  boundary  layer.  These 
integrations  yield 


(3.14) 


4 

i 


for  the  heat  flux  distribution.  Since  the  fluid  dynamics  data  were  taken 
under  isothermal  conditions,  Equation  (3.14)  assumes  p = p^Cx) , while 
Equation  (3.15)  retains  p = p(x,y).  (Although  these  variations  were 
included  in  :he  analysis,  numerically  they  were  insignificant  in  all  cases.) 

The  shear  stress  and  heat  flux  contain  both  turbulent  and  laminar 
contributions  and  may  be  written  as 


T = -p  u'v'  + p V 
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(3.16) 


and 


q"  - PCpv ' t ' - fcfi 


(3.17) 


If  the  turbulent  contributions  are  modeled  using  eddy  diffusivities  for 
momentum  and  heat,  Equations  (3.16)  and  (3.17)  become 
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and 


• AT 

q"  = - P Cp  (eH  + a)  ^ 
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where 


- u'v'  = e 
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and 
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Alternatively,  if  the  turbulent  shear  stress  is  modeled  using  the 
mixing  length  approach,  a mixing  length  8.  may  be  defined  as 


z 


m 


au 


8y 


(3.22) 


or,  using  Equation  (3.20) 


(3.23) 


The  turbulent  Prandtl  number  is  defined  as  the  ratio  of  the  eddy 
diffusivities  for  momentum  and  heat: 

PrT  = em/eH  (3.24) 

In  order  to  demonstrate  the  consistency  of  the  Pr-p  results  calculated  using 
the  method  outlined  above  with  the  Pr^,  data  obtained  by  Pimenta  [2]  from 
measurements  of  u'v' , v't',  and  dT/dU,  values  of  Pr^.  were  calculated  using 
the  present  method  for  the  unaccelerated,  unblown  = 89  ft/sec  case 
reported  by  Pimenta.  Results  of  this  calculation  are  compared  in  Figure  3-38 
with  the  measured  values  of  PrT  reported  by  Pimenta.  The  two  methods  give 
results  which  agree  well  in  the  inner  region.  The  calculated  data  are  very 
uncertain  in  the  outer  region  where  ^ and  ^ approach  zero  since  the 
uncertainty  in  the  numerical  calculations  of  the  derivatives  approaches 
infinity  as  y 6.  Pimenta  avoided  this  increase  in  uncertainty  by  calcu- 
lating dT/dU  from  the  linear  (Tw-T)/ (T^-T^)  vs.  U/U^  plots  discussed 
previously. 
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Calculated  values  of  Pr,p  for  the  four  acceleration  cases  of  this  study 
are  presented  in  Figure  3-39.  Also  shown  are  the  bounds  on  the  smooth  wall 
acceleration  data  reported  by  Kearney  for  K < 2.5  x 10  6 and  the  calculated 
data  for  Pimenta's  Kr  = 0 case.  The  rough  wall  data  lie  at  the  lower  edge 
of  the  smooth  wall  data  range.  It  appears  from  the  present  data  that  the 
use  of  a constant  Pr^  = 0.7  - 0.8  would  be  a reasonable  assumption  in  a 
prediction  method  modeling  accelerated  flow  over  the  present  rough  surface. 
Comparison  of  the  present  data  with  those  of  Pimenta  indicates  that  for  fully 
rough  flow  the  turbulent  Prandtl  number  is  decreased  slightly  by  acceleration. 

Mixing  lengths  calculated  using  Equation  (3.23)  are  presented  in  Figure 



3-40.  In  the  determination  of  l,  values  of  u'v'  calculated  from  Equations 
(3.14)  and  (3.16)  were  used  since  comparison  with  the  measured  values  showed 
agreement  within  a few  percent  and  the  calculations  yielded  values  of  u'v' 
closer  to  the  surface  than  were  possible  to  obtain  with  the  probes.  The 
plot  shows  that  in  the  inner  region  £/(y+Ay)  is  slightly  lower  than  the 
unaccelerated  values  of  0.40  - 0.41  found  by  Pimenta.  The  behavior  of  Jt/6 
in  the  outer  region  is  in  agreement  with  that  observed  by  Pimenta  for  Kr  = 0. 

Calculations  of  the  nondimens ional  shear  stress,  x + , and  nondimens ional 
heat  flux,  Q+,  are  presented  in  Figure  3-41  for  the  Kr  = 0.29  x 10  F = 0 
and  0.0039  runs.  The  trends  observed  are  similar  to  those  for  smooth  wall 
accelerated  flows  [25],  indicating  that  the  effect  of  roughness  on  these 
nondimens ional  distributions  is  small.  The  turbulent  contributions  to  x+ 
and  Q+  are  also  shown  (denoted  with  the  subscript  T) . Comparison  of  the 
total  values  with  the  turbulent  ones  shows  that  the  laminar  contribution  to 
both  shear  stress  and  heat  flux  was  extremely  small  throughout  the  region  of 
measurement . 
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The  expression  for  q”  is  the  one  normally  found  by  assuming  all  terms  to  be 
zero.  However,  in  Equation  (3.26)  for  tc  only  the  first  two  terms  on  the 
right  hand  side  are  the  ones  normally  retained  for  the  Couette  flow  approxi- 
mation. The  additional  term  can  be  viewed  as  a correction  term  for  the 
effects  of  the  pressure  gradient.  Thus  the  expression  for  tc  might  be  more 
accurately  labeled  a "near-wall"  approximation  rather  than  a Couette  flow 
approximation.  In  any  case,  it  can  be  observed  from  the  figure  that 
Equations  (3.25)  and  (3.26)  provide  accurate  representations  of  the 
behavior  of  t+  and  Q+  in  the  region  very  near  the  wall. 

Values  of  v' t'  for  the  two  unblown  equilibrium  runs  calculated  from 
Equations  (3.15)  and  (3.17)  are  shown  in  Figure  3-42  as  v't' /UtTtvs  y/6. 
Comparison  of  these  values  with  the  unaccelerated  data  measured  by  Pimenta  [2] 
and  Orlando  [34]  on  rough  and  smooth  walls,  respectively,  indicates  that  the 
distribution  of  v't* AJtTt  is  independent  of  acceleration  and  surface  condi- 
tion, at  least  within  the  ranges  of  the  data  available. 

Figure  3-43  shows  results  of  calculation  of  the  turbulent  kinetic  energy 
production  for  Pimenta's  [2]  zero  pressure  gradient  data  and  the  present 
Kr  = 0.15  x 10  3 data.  From  Equation  (3.11),  the  production  term  can  be 
written  (using  the  standard  boundary  layer  assumptions)  as 


i 
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(3.27) 
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P = -U'V'  ~ U'  7T-  + v' 
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The  second  and  third  terms  are  normally  neglected  in  zero  pressure  gradient 

flows.  In  the  calculations  presented  in  Figure  3-43,  the  last  term  in 

2 

Equation  (3.27)  was  neglected  since  measurements  of  u'  were  made  much 

2 

closer  to  the  wall  than  those  of  v'  . Thus  the  results  shown  present  an 

91) 

upper  bound  on  the  effect  of  the  pressure  gradient  through  the  ^ terms. 

In  a boundary  layer  subjected  to  a favorable  pressure  gradient,  the 


second  term  in  (3.27)  is  negative  and  thus  appears  as  a sink  for  turbulent 

2 

kinetic  energy.  Hinze  [22]  notes  that  one  should  expect  a decrease  in  q 

as  a result.  Such  a decrease  in  turbulence  energy  was  noted  in  the  present 

accelerated  data  (see  Figure  3-31,  for  example).  However,  the  results  shown 

in  Figure  3-43  indicate  that  for  the  present  data  the  production  is  decreased 

with  acceleration  primarily  because  of  changes  in  the  distributions  of 
9U 


-u'v'  and  , while  the  sink  term  remains  of  negligible  magnitude, 
ay 


3.3  Heat  Transfer  Predictions  and  Supplementary  St  Data 

The  Stanton  number  data  discussed  in  previous  sections  were  all  for 

constant  F,  constant  Tw  boundary  conditions.  In  the  period  following  the 

investigation  reported  by  Pimenta  [2]  and  in  the  course  of  the  present  study, 

additional  St  cases  were  run  for  conditions  of  variable  of  T and  F.  A 

w 

prediction  method  for  rough  wall  heat  transfer  with  variable  velocity,  wall 
temperature  and  blowing  was  developed  using  these  data.  In  this  section, 
the  prediction  method  is  described  and  some  typical  comparisons  of  data  and 
predictions  are  presented.  Additional  cases  of  St  with  variable  boundary 
conditions,  not  presented  in  the  Figures,  are  tabulated  in  Appendix  IV. 
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For  smooth  walls  the  variable  wall  temperature  case  was  dealt  with  by 


Reynolds,  et  al.  139],  using  superposition  based  on  a kernel  function 
describing  the  downstream  effects  of  a step  change  in  wall  temperature. 

Whitten  [40]  extended  this  to  include  variable  blowing  and  Orlando  [34]  used 
this  same  method  for  variable  wall  temperatures  in  adverse  pressure  gradients. 
However,  Orlando  found  that  smooth  wall  Stanton  numbers  were  not  affected  by 
the  adverse  pressure  gradients  investigated  when  presented  in  enthalpy 
thickness  coordinates. 

It  was  shown  by  Healzer  [1]  and  Pimenta  [2]  over  a range  of  free 
stream  velocities  that  for  a fully  rough  turbulent  boundary  layer  flow  on 
the  present  surface,  with  constant  U^,  Tw  and  F 


from  which  it  follows  that 


where  r is  the  radius  of  the  spheres  comprising  the  test  surface.  Thus  there 
is  a unique  curve  of  St  vs  x/r  for  each  value  of  F,  so  long  as  both  Uot  and 
Tw  are  uniform.  The  present  data  for  constant  and  AT  confirm  this,  being 
well  represented  by  the  interpolation  expression 


where 


A = -1.36  + 48.2  F 


B = -0.61  - 57.4  F 
C = 0.0675  + 3.69  F 


for  0 < F < 0.004.  Figure  3-44  shows  the  StQ  data  and  interpolation  curves 
from  Equation  (3.30)  for  three  values  of  F.  The  F = 0.002  data  are  from 
Pimenta's  [2]  study  and  were  untripped,  while  the  data  for  F = 0 and  0.0039 
are  from  the  present  study  which  did  use  a boundary  layer  trip. 

Now  consider  the  case  of  an  unheated  starting  length  with  constant  F 
and  UM.  Reynolds,  et  al.  [39],  showed  that  for  turbulent  flow  over  a 
smooth  wall  with  F = 0 the  Stanton  number  downstream  of  a step  increase  in 
wall  temperature  could  be  predicted  by 


St 

Sto 


x > £ 


(3.31) 


where  £ is  the  x-position  of  the  temperature  step,  m = 9/10  and  n = -1/9. 
This  expression  was  then  used  as  the  kernel  function  in  a superposition 
integral  for  predictions  with  arbitrarily  varying  AT,  as  follows: 


d[AT(S)l  , 

— a? — ^ 


(3.32) 


A kernel  function  of  the  same  type  can  be  developed  for  the  present 
rough  wall  unheated  starting  length  results,  and  the  data  are  well  represented 
using  m = 1 and  n = -0.22.  Thus,  for  fully  rough  turbulent  flow  with  unheated 
starting  length  and  U^,  F and  AT  constant 


(If'22 
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x > £ 


(3.33) 


where  St  is  evaluated  at  the  same  F and  x as  St  . Figure  3-45  shows  the 
unheated  starting  length  St  data  and  curves  evaluated  using  Equation  (3-33) 
for  two  cases  with  F = 0 and  one  case  with  F = 0.0039. 

As  shown  previously,  for  a fully  rough  turbulent  boundary  layer  a 
positive  value  of  dll^/dx  gives  a higher  Stanton  number  than  the  constant 
velocity  case  at  the  same  x or  A?.  Thus,  the  response  of  Stanton  number  to 
positive  dl) ^dx  is  the  same  as  to  positive  values  of  dT ^dx.  Also,  consider- 
ation of  the  energy  integral  equation  (with  constant  properties)  in  the  form 

st  4 F * or  a?  (42U.4T)  <3-31> 

00 

shows  that  the  variables  and  AT  always  appear  in  the  product  form  as 
(U^AT).  (Of  course,  Equation  3.34  is  a conservation  equation  obeyed  by  both 
smooth  and  rough  wall  flows.  The  different  response  of  smooth  and  fully 
rough  layers  to  positive  dU J dx  is  evidently  related  to  the  effect  of  the 
viscous  sublayer  present  in  smooth  wall  layers  on  the  rate  equations.) 

These  observations  lead  to  the  following  proposed  prediction  method  for 
cases  of  variable  U^,  AT  and  F in  a fully  rough  flow: 

1~~22  dT  KW  *1(01  x>(  (3.35) 

(1-5/x) 

where  St  is  evaluated  at  the  same  x and  F as  St  . The  predictions  which 
follow  were  obtained  by  numerically  integrating  Equation  (3.35).  For  these 
calculations,  the  integral  was  expanded  to  a sum  of  two  integrals  and  the 
dUm(0  term  was  approximated  by  assuming  a linear  variation  of  with  £ 
across  each  integration  step. 

It  was  found  that  Equation  (3.35)  worked  well  except  in  cases  where 

there  were  steps  in  F.  Obviously,  St  cannot  change  instantaneously  after  a 
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step;  therefore,  modification  of  the  method  is  required  to  account  for  the 


"lag"  in  St  after  such  a step.  A simple  and  satisfactory  approach  is  to 
define  a new  StQ  which  is  modified  by  the  kernel  function,  giving 

N 

St*  = St  - ASt 
O O o 

i=l 

where  StQ.  = StQ(8*)  - Sto(iL)  and  8-  is  the  position  of  the  ith  step.  Thus, 
when  there  are  steps  in  F,  use  of  St*  in  place  of  StQ  in  Equation  (3.35) 
accounts  for  the  "lag”  in  St  caused  by  the  step. 

Figure  3-46  presents  data  and  predictions  for  constant  Uro,  F = 0 and  a 
bilinear  variation  of  wall  temperature;  Figure  3-47  shows  cases  for  the 
K = .28  x 10  ^ nonequilibrium  run  both  with  and  without  unheated  starting 
length;  and  in  Figure  3-48  data  and  predictions  for  an  unblown,  Kr  = .50  x 10 
equilibrium  run  are  presented.  In  all  three  figures  the  agreement  between 
the  data  and  predictions  is  very  good. 

Figures  3-49  and  3-50  present  data  from  two  cases  designed  to  provide  a 
test  of  any  prediction  method  proposing  to  calculate  heat  transfer  in  the 
fully  rough  turbulent  boundary  layer.  The  flow  is  accelerated  arbitrarily 
and  subjected  to  a step  in  F,  followed  by  a variable  F and  then  a step  back 
to  F = 0.  Figure  3-49  presents  the  AT  constant  case,  while  Figure  3-50 
shows  results  when  a step  in  AT  is  imposed  in  the  region  of  variable  F.  In 
both  figures  the  dashed  line  is  the  prediction  using  Equation  (3.35),  while 
the  solid  line  shows  the  prediction  including  the  modification  for  the  steps 
in  F (Equation  (3.36)).  It  is  evident  that  the  lag  introduced  by  the  steps 
in  F must  be  taken  into  account,  but  once  this  is  done  the  present  predic- 
tions are  in  excellent  agreement  with  the  data. 
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Figure  3-1.  Zero  Pressure  Gradient  Stanton  Number  Data 


Figure  3-2.  Zero  Pressure  Gradient  Skin  Friction  Coefficient  Data 
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Figure  3-3.  Typical  Mean  Velocity  Profile  for  Fully 
Rough,  Zero  Pressure  Gradient  Flow 


Figure  3-4.  Conponents  of  Turbulent  Kinetic  Energy  for 
Fully  Rough,  Zero  Pressure  Gradient  Flow 
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Figure  3-5.  Suninary  Data  for  Kr  = 0.15  x 10 

F = 0 Equilibrium  Acceleration  Case 


51 


0 10  20  30  40  50  60  70 

x (in.) 


Figure  3-8.  Summary  Data  for  K = 0.28  x 10  , 

F = 0 Nonequilibrium  Acceleration  Case 
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Figure  3-9.  Comparison  of  Enthalpy  Thickness  Data  and  Solution 
for  Equilibrium  Enthalpy  Thickness  Behavior 
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Figure  3-10.  Equilibrium  Acceleration  Stanton  Number  Data 
vs  (Enthalpy  Thickness/Sphere  Radius) 
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Roughness  Reynolds  Number  Variation  with  Axial  Distance 
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Figure  3-18.  Mean  Velocity  Profiles  Illustrating  Lack  of 

Similarity  in  Flow  Direction  (K  = 0.28  x 10'°) 
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Figure  3-20.  Mean  Velocity  Profile  Plotted  Using  Shifted 
Wall  Position  (Kr  = 0.29  x 10" F = 0) 
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3-21.  Comparison  of  Blown  and  Unblown  Mean  Velocity  ^ 
Profiles  at  Same  Axial  Position  (K  = 0.29  x 10  ) 


3-22.  Nondimens ional  Mean  Temperature  Profiles  Illustrating 
Similarity  in  Flow  Direction  (K  = 0.15  x 10'^) 
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Figure  3-23.  Nondimens ional  Mean  Temperature  Profiles  Illustrating 
Similarity  in  Flow  Direction  (Kr  = 0.29  x 10'3,  F = 0) 
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Figure  3-24.  Nondimens ional  Mean  Temperature  Profiles  Illustrating 

Similarity  in  Flow  Direction  (K  = 0.29  x 10' 3,  F = 0.00 
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Nondimens ional  Mean  Temperature  Profiles  for 
Nonequilibrium  Acceleration  Run  (K  = 0.28  x 10 


Figure  3-25 
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Figure  3-26.  Nondimens ional  Mean  Temperature  Plotted  versus 
Nondimens ional  Mean  Velocity  (K  = 0.15  x 10*3 
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-27.  Comparison  of  Blown  and  Unblown  Nondimeioional  Mean  ^ 

Temperature  Profiles  at  Same  Axial  Position  (Kr  = 0.29  x 10  ) 
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’ irbulcnt  Kinetic  Fnergy  Components  Illustrating 
low  Direction  (Kr  = 0.15  x 10" 3) 
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Figure  3-29.  Profiles  of  Reynolds  Shear  Stress  Compared 
at  TV/o  Axial  Positions  (K^  = 0.15  x 10" 3) 
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Figure  3-30.  uf  /U  Profiles  Compared  for 
K = 0 and  K = 0.15  x 10' 3 


Figure  3-31.  Profiles  of  Components  of  Turbulent  Kinetic 
Energy  Compared  for  Kf  = 0 and  K^  = 0.15  x 10 


Figure  3-32.  Profiles  of  Components  of  Turbulent  Kinetic 

Energy  Compared  for  K = 0 and  K = 0.29  x 10“’ 
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Figure  3-33.  Profiles  of  Components  of  Turbulent  Kinetic  Energy 
Compared  for  Kr  = 0 and  K = 0.28  x 10‘6 
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Figure  3-34.  Profiles  of  Components  of  Turbulent  Kinetic  Energy 

Compared  for  K = 0,  F = 0.0039  and  K-  0.29  x 10‘3 
F = 0.0039  r r 
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Figure  3-35.  Comparison  of  Blown  and  Unblown  Profiles 
of  u'^/U^  for  Unaccelerated  (Kr  = 0) 
and  Accelerated  (Kr  = 0.29  x 10  Cases 


Figure  3-36.  Comparison  of  Reynolds  Shear  Stress 

Profiles  for  K_  = 0 and  K_  = 0.15  x 10 
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Figure  3-38.  Comparison  of  Turbulent  Prandtl  Numbers  Measured  by 

Pimenta^  and  Calculated  by  Present  Method  for  Kr  = 0 


1 i I TTTTT]  i — i — I I I II  1 1 i — i — I II  II  i 

Kr  = 0 F = 0 (PIMENTA) 

• MEASURED  BY  PIMENTA 
o CALCULATED  BY  PRESENT  METHOD 


00°°°0oo0°°Oo S8S*°°8J.  . . 

°o0  •- 

o 

o 

o 

o- 

J 1 — LI  LLLlJ 1 I — I I I I I I I I I I II  II  I 


2.4 

2.0 

1.6 

PrT 

1.2 

.8 


.4 
0 

.01  .10  1.0 

y/8 

Figure  3-39.  Turbulent  Prandtl  Number  Distributions  for  the 
Acceleration  Cases  of  the  Present  Study 
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Figure  3-40.  Mixing  Length  Distributions  for  the  Equilibrium 
Acceleration  Cases  of  the  Present  Study 
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Figure  3-41.  r and  Q Distributions  for  the  K = 0.29  x 10  , F = 0 and  F = 0.0039  Runs 


Figure  3-42.  Comparison  of  the  Unblown  v't'/UTTT  Profiles 

for  the  Equilibrium  Accelerations  of  the  Present 
Study  with  the  Unaccelerated  Profiles  for  Smooth 
and  Rough  Wall  Layers 
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Figure  3-43.  Calculated  Distribution  of  the  Nondimens ional 

Turbulent  Kinetic  Energy  Production  for  K = 

and  K = 0.15  x 10'3 
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Figure  3-44.  Comparison  of  St  Data  and  Interpolation 
Expression  for  U^,  Tw  and  F Constant 
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Figure  3-45.  Comparison  of  Unheated  Starting  Length  St  Data 

with  Results  of  Kernel  Solution  for  U and  F Constant 


Comparison  of  St  Data  and  Pre-  Figure  3-47.  Comparison  of  St  Data  and  Pre- 
diction for  a Bilinear  Variation  dictions  for  K = 0.28  x 10'6, 

with  Uro  Constant,  F = 0 F = 0 Run  Both  with  and  Without 

an  Unheated  Starting  Length 
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(Figure  3-48.  Comparison  of  St  Data  and  Prediction  for 

Kr  = 0.50  x 10" 3,  F = 0 Run  with  Tw  Constant 
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CHAPTER  4 


EFFECTS  OF  ACCELERATION  ON  THE  FULLY  ROUGH  TURBULENT  BOUNDARY  LAYER 


In  the  previous  chapter  the  various  sets  of  data  taken  in  this  study 
were  presented  and  some  details  of  the  boundary  layer  behavior  indicated  by 
the  data  were  discussed.  A more  general,  integrated  discussion  of  the 
effects  of  acceleration  on  the  fully  rough  turbulent  boundary  layer  is 
presented  in  this  chapter. 

In  smooth  wall  turbulent  boundary  layers,  a viscous  sublayer  develops 
which,  in  a sense,  provides  a "buffer"  between  the  surface  and  the  more 
energetic  turbulent  fluid  in  the  outer  portion  of  the  layer.  For  flow  on  a 
rough  surface,  once  the  boundary  layer  is  in  the  fully  rough  state  the 
viscous  sublayer  has  been  effectively  destroyed,  at  least  above  the  crests 
of  the  roughness  elements.  The  outer  flow,  therefore,  interacts  with  near- 
wall  fluid  of  higher  momentum,  energy  and  vorticity  than  in  the  smooth  wall 
case.  Thus,  the  fully  rough  flow  near  the  wall  is  characterized  by  shorter 
time  scales  and  larger  velocity  scales  than  for  a smooth  wall  flow  and  is 
unaffected  by  molecular  viscosity  down  to  the  tops  of  the  roughness  elements. 
Among  the  elements,  the  pressure  fluctuations  are  thought  to  be  more  important 
than  the  viscous  forces  [2].  Blake  [21]  found  that  pressure  fluctuation 
intensity  near  the  wall  increased  with  roughness  and  scaled  with  the  wall 
shear  stress  when  compared  to  smooth  wall  values. 

In  view  of  these  characteristics,  it  is  not  surprising  that  the  proper 


acceleration 


/ r dM 

ion  parameter  for  fully  rough  flow  lKr  = ^ ) is  independent 
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of  molecular  viscosity  and  contains  a length  scale  dependent  on  the  wall 
roughness  elements.  As  shown  in  Chapter  3,  imposition  of  a Kj,  = constant 
acceleration  on  a fully  rough  flow  with  F constant  and  > 0 results  in  a 
boundary  layer  which  approaches  an  asymptotic  condition  for  which  C^/2,  62, 
H,  0,  and  G are  all  constant  with  x.  This  type  of  boundary  layer  satisfies 
the  conditions  found  by  Rotta  [16]  for  an  exactly  equilibrium  flow,  i.e.. 


and 


C^/2  = constant 


d6j/dx  = constant  (=0) 


= constant 


In  the  equilibrium  accelerated  fully  rough  layer,  the  mean  velocity 
profiles  become  similar,  and  the  profiles  of  the  Reynolds  stress  tensor 
components  also  approach  a similar  behavior.  The  roughness  Reynolds  number, 
Re^,  increases  throughout  the  acceleration  region,  indicating  an  evolution 
toward  a "rougher"  layer  rather  than  a transitionally  rough  or  smooth 
behavior.  This  is  opposite  to  the  behavior  of  the  accelerated  smooth  wall 
turbulent  boundary  layer,  which  evolves  toward  a state  resembling  laminar 
flow.  These  differences  in  behavior  are  consistent  with  the  Stanton  number 
trends  observed  with  acceleration.  In  the  fully  rough  layer,  Stanton  numbers 
are  increased  with  acceleration  compared  to  unaccelerated  values  at  the  same 
x,  ^2»  or  Re£  • In  the  smooth  wall  layer,  however,  Stanton  numbers  are  either 
unchanged  or  are  decreased  (depending  on  the  strength  of  the  acceleration) 
compared  with  unaccelerated  values  at  the  same  position  or  Reynolds  number. 


.1 


j 


0 j 
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The  virtual  position  of  the  wall,  found  using  the  method  suggested  by 
Monin  and  Yaglom  [30]  and  based  on  the  mean  velocity  profiles,  is  unaffected 
by  either  acceleration  or  axial  position  x.  The  wall  shift,  Ay,  reported  by 
Pimenta  [2]  for  the  present  surface  with  Kr  = 0 remained  valid  under  all 
conditions  investigated  in  this  study. 

Although  the  accelerated  (T  -Tj/H^-Tj  vs  U/U^  data  do  not  exhibit  the 
linearity  observed  in  = 0 data  [ 2 ] , a temperature  "jump"  at  the  extra- 
polated U/U^  = 0 point  is  still  indicated.  Therefore,  there  are  different 
apparent  wall  positions  for  the  mean  velocity  and  mean  temperature  profiles 
in  both  zero  pressure  gradient  and  accelerated  fully  rough  flows,  with  the 
virtual  wall  indicated  by  the  mean  temperature  profiles  being  further  below 
the  crests  of  the  roughness  elements  than  the  apparent  wall  determined  from 
the  mean  velocity  profiles. 

The  turbulence  field  is  dramatically  affected  when  a fully  rough 

turbulent  boundary  layer  is  accelerated.  The  nondimensionalized  turbulent 
2 2 2 2 

kinetic  energy  (q  /U^  or  q /UT)  is  substantially  reduced  over  the  entire 
layer  compared  to  the  unaccelerated  case.  This  is  the  same  trend  that  can 
be  inferred  from  the  data  for  accelerated  smooth  wall  layers.  Although  no 


measurements  of  v'  and  w'  in  accelerated  smooth  wall  layers  are  known  to 
— 1 2 

the  author,  u'  /\Jm  profiles  are  decreased  by  acceleration  [5],  and  a decrease 

in  bursting  has  been  noted  with  acceleration  [36].  The  present  data  also 

show  that  an  accelerated  fully  rough  layer  is  much  more  anisotropic  in  the 

inner  region  than  is  a K = 0 fully  rough  layer.  This  is  probably  due  to 

r / 

the  influence  of  acceleration  on  the  pressure  field-velocity  field  interaction 

terms  Ip'  , p'  1 which  transfer  energy  from  the  u'  to  the  v'  and  w' 

' 2 
components  of  q . 
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In  the  accelerated  layer  with  blowing  investigated  in  this  study,  the 
2 2 

decrease  in  q /Um  in  the  outer  region  observed  in  the  unblown  cases  was 
noted,  but  the  increased  anisotropy  in  the  inner  region  observed  for  F = 0 
was  not  seen  for  F = 0.0039.  This  trend  was  also  noted  in  the  Kr  = 0 fully 
rough  data  of  Pimenta  [2],  who  found  that  blowing  produced  a more  isotropic 
turbulence  field  than  that  of  an  unblown  layer.  It  thus  appears  that,  in 
the  present  data,  the  effects  of  acceleration  and  blowing  on  the  turbulence 
field  in  the  inner  region  of  the  layer  are  approximately  equal  and  opposite. 

The  values  of  the  Reynolds  shear  stress  correlation  coefficients  found 

2 

in  this  study  are  in  agreement  with  the  values  (Ruv  =»  0.45,  R^  = 0.15) 
previously  reported  for  smooth  wall  [32,33,34]  and  = 0 rough  wall  [2] 
turbulent  boundary  layers.  This  observation  indicates  a universal  mechanism 
in  the  bursting  and  decay  process.  Although  the  magnitudes  of  the  components 
of  the  Reynolds  stress  tensor  are  dependent  on  surface  condition,  blowing 
and  pressure  gradient,  the  constancy  of  the  correlation  coefficients 
indicates  a universal  and  constant  relationship  among  these  components. 

There  are  indications  from  the  present  data  that  in  a fully  rough 
turbulent  boundary  layer  with  Kr,  F,  and  T constant,  the  layer  approaches 
an  equilibrium  state  in  the  thermal  sense  also.  Unfortunately,  the  range  of 
the  present  thermal  data  is  not  large  enough  to  allow  a definitive  conclusion 
in  this  regard.  It  is  proposed  that  in  such  an  equilibrium  state  Stanton 
number  would  become  constant,  enthalpy  thickness  would  approach  an  asymptote, 
and  nondimensional  temperature  profiles  would  be  similar  in  the  flow 
direction.  In  the  present  data  the  fluid  dynamic  field  (which  is  a major 
influence  on  the  development  of  the  thermal  field)  is  in  an  equilibrium 
state  described  previously,  and  the  mean  temperature  profiles  are  similar  in 
the  flow  direction.  However,  due  to  the  relatively  small  range  of 
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covered  in  the  accelerating  regions,  it  is  not  possible  to  show  that  the 
present  Stanton  number  data  are  constant  rather  than  weak  functions  of 
enthalpy  thickness. 

As  in  the  case  of  the  Reynolds  shear  stress  correlation  coefficients, 
the  distribution  of  the  turbulent  heat  flux  correlation  coefficient 
v't'  / UtTt  calculated  from  the  data  of  the  present  study  is  identical 
with  the  profiles  reported  for  zero  pressure  gradient  rough  wall  [2]  and 
smooth  wall  layers  [34]. 

Results  of  turbulent  Prandtl  number  calculations  based  on  data  of  the 
present  experiment  indicate  an  approximately  constant  value  of  0.7  - 0.8 
across  the  layer.  This  value  is  lower  than  turbulent  Prandtl  numbers 
reported  for  the  unaccelerated  fully  rough  case  [2],  which  vary  from  ~ 1.0 
in  the  near  wall  region  to  ~ 0.7  - 0.8  in  the  outer  region.  The  present 
values  of  Pr^,  are  in  the  lower  range  of  the  data  reported  by  Kearney  [25] 
for  smooth  wall  accelerated  layers. 

As  shown  in  Section  3.4,  the  response  of  Stanton  number  in  fully  rough 
flow  to  either  acceleration  or  variable  wall  temperature  can  be  represented 
using  the  same  kernel  function,  |l  - 
the  absence  of  a viscous  sublayer,  which  leads  to  greater  turbulent  mixing 
and  more  vigorous  interaction  between  the  wall  region  and  outer  flow  [38.2]. 
By  contrast,  in  a smooth  wall  layer  the  viscous  sublayer  damps  or  buffers 
the  interaction  between  the  wall  and  the  fully  turbulent  layer,  and  the 
response  of  Stanton  number  is  different  for  variations  in  velocity  and  wall 
temperature . 

The  fact  that  variations  in  velocity  and  temperature  have  equivalent 
effects  on  Stanton  number  and  that  the  process  can  be  modeled  in  the  manner 
described  above  indicates  that  turbulent  Prandtl  number  should  be  unity. 


*) 


-u.zz 


This  behavior  is  related  to 
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This  information  complements  the  results  of  the  PrT  calculations  described 
above,  since  the  calculated  values  are  of  order  one. 
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CHAPTER  5.  CONCLUSIONS 


All  statements  in  this  chapter,  unless  otherwise  indicated,  refer  to  the 
flow  of  a fully  rough  turbulent  boundary  layer  over  a three-dimensional, 
densely  packed,  uniformly  rough  test  surface  such  as  that  used  in  this  study. 
Based  on  the  discussion  in  the  previous  chapters,  the  important  results  and 
conclusions  of  the  study  are: 

1.  The  proper  acceleration  parameter  for  use  with  fully  rough  flow 

dU 

T*  OO 

is  Kr  = o — gj-  » where  r is  a characteristic  roughness  length. 

OO 

2.  In  a constant  Kr  acceleration  with  F > 0 and  constant,  the  fully 
rough  layer  develops  toward  an  equilibrium  state  where  C^/2,  62,  H, 

8 and  G are  all  constant.  Both  mean  velocity  profiles  and  the 
components  of  the  Reynolds  stress  tensor  approach  similarity  in  the 
flow  direction. 

3.  Although  inherent  uncertainty  in  the  Stanton  number  data  and  the 

restricted  length  of  the  acceleration  region  prevent  a definite 

conclusion,  the  present  thermal  data  indicate  the  possibility  that 

a layer  with  K , T and  F(>  0)  constant  approaches  an  equilibrium 
I w 

state  in  the  thermal  sense,  also.  Such  a state  would  be  character- 
ized by  Stanton  number  becoming  constant,  enthalpy  thickness 
approaching  an  asymptote,  and  temperature  profiles  being  similar  in 


the  flow  direction. 
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4.  Stanton  numbers  increase  with  acceleration  compared  to  zero  pressure 
gradient  values  at  the  same  position  or  enthalpy  thickness.  This 
behavior  is  quite  different  from  that  of  a smooth  wall  layer,  where 
Stanton  numbers  are  unchanged  for  small  values  of  K and  decrease 
with  increasing  K compared  to  zero  pressure  gradient  data  at  the 
same  position  or  Reynolds  number. 

5.  Roughness  Reynolds  number  increases  in  a region  of  acceleration, 
indicating  that  the  fully  rough  layer  does  not  tend  toward  the 
transitionally  rough  or  smooth  wall  state  when  accelerated. 

6.  For  F = 0,  acceleration  decreases  the  turbulent  kinetic  energy 
throughout  the  boundary  layer,  and  in  the  inner  region  the  turbulence 
is  much  more  anisotropic  than  in  the  Kr  = 0 layer. 

7.  In  the  blown  layer,  acceleration  decreases  the  turbulent  kinetic 

energy  in  the  outer  region  of  the  layer  and  substantially  alters  the 

2 2 Z 

shape  of  the  profiles  of  u'  , v'  , w'  compared  with  the  blown, 
unaccelerated  boundary  layer  (Figure  3-34). 

8.  The  values  of  the  Reynolds  shear  stress  coefficients  obtained  in 
this  study  are  the  same  as  those  reported  for  smooth  wall  flows 

and  Kr  = 0 fully  rough  flows,  indicating  that  these  values  are  truly 
universal. 

9.  The  profiles  of  v’t’/U  T^  calculated  from  the  present  accelerated 
data  are  in  good  agreement  with  those  reported  for  smooth  and  rough 
wall  zero  pressure  gradient  layers.  It  thus  appears  that  the 
turbulent  heat  flux  profile,  nondimensionalized  in  this  manner,  is 
independent  of  surface  condition  and  favorable  pressure  gradient. 

10.  In  a fully  rough  flow,  the  response  of  Stanton  number  to  either 

variable  wall  temperature  or  acceleration  can  be  represented  using 

/ g\-0.22 

the  same  kernel  function,  ll  - -) 
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11.  From  (10)  above  one  would  expect  a turbulent  Prandtl  number  of 
approximately  unity.  The  values  of  PrT  calculated  from  the  present 
data  support  this,  having  an  approximately  constant  value  of  0.7  - 
0.8  across  the  boundary  layer. 

12.  The  virtual  origin  of  the  wall  (based  on  mean  velocity  profiles)  is 
independent  of  acceleration  and  position  in  the  flow  direction. 


REFERENCES 


1.  Healzer,  J.  M. , Moffat,  R.  J.,  and  Kays,  W.  M. , ’’The  Turbulent  Boundary 

Layer  on  a Rough,  Porous  Plate:  Experimental  Heat  Transfer  with  Uniform 

Blowing,"  Report  HMT-18,  Thermosciences  Division,  Dept,  of  Mech.  Eng., 
Stanford  University,  1974. 

2.  Pimenta,  M.  M.,  Moffat,  R.  J.,  <*nd  Kays,  W.  M. , "The  Turbulent  Boundary 

layer : An  Experimental  Study  of  the  Transport  of  Momentum  and  Heat  with 

the  Effect  of  Roughness,"  Report  HNfT-21,  Thermosciences  Division,  Dept, 
of  Mech.  Eng.,  Stanford  University,  1975. 

3.  Nikuradse,  J. , "Stromugsgestze  in  rauhen  Rohren,"  VDI  Forschungsheft, 

No.  361,  1950,  English  Translation,  NACA  TM  1292. 

4.  Kays,  W.  M.,  and  Moffat,  R.  J.,  "The  Behavior  of  Transpired  Turbulent 
Boundary  Layers,"  Report  U4T-20,  Thermosciences  Division,  Dept,  of 
Mech.  Eng.,  Stanford  University,  1975. 

5.  Badri  Narayanan,  M.  A.,  and  Ramjee,  V.,  "On  the  Criteria  for  Reverse 
Transition  in  a Two-Dimensional  Boundary  Layer  Flow,"  Jn.  Fluid  Mech., 

Vol.  35,  pp.  225-241,  1969. 

6.  Reshotko,  M.  Boldman,  D.  R. , and  Ehlers,  R.  C.,  "Heat  Transfer  in  a 60° 
Half-Angle  of  Convergence  Nozzle  with  Various  Degrees  of  Roughness," 

NASA  TN  D-5887 , 1970. 

7.  Banerian,  G. , and  McKillop,  A.  A.,  "The  Effects  of  Surface  Roughness  in 
Nozzles  on  Heat  Transfer,"  Proc.  of  the  Fifth  International  Heat  Transfer 
Conference , V.  II,  pp.  234-238,  1974. 

8.  Chen,  K.  K. , "Compressible  Turbulent  Boundary- Layer  Heat  Transfer  to 
Rough  Surfaces  in  Pressure  Gradient,"  AIAA  J. , 10,  pp.  623-629,  1972. 

9.  Blackwell,  B.  F. , Kays,  W.  M. , and  Moffat,  R.  J.,  "The  Turbulent 

Boundary  Layer  on  a Porous  Plate:  An  Experimental  Study  of  the  Heat 

Transfer  Behavior  with  Adverse  Pressure  Gradients,"  Report  No.  HMT-16, 
Thermosciences  Division,  Dept,  of  Mech.  Eng.,  Stanford  University,  1972. 

10.  Rotta,  J.,  "Uber  die  Theorie  der  turbulenten  Grenzschichten," 

Mitteilungen  aus  dem  Max-Planck-Institut  fur  Stromungsforschung 
(Gottingen),  No.  1,  1950.  English  translation  NACA  TM  1344,  1953. 

11.  Clauser,  F.  H. , "Turbulent  Boundary  Layers  in  Adverse  Pressure  Gradients," 
Jn.  Aero.  Sci.,  Vol.  21,  pp.  91-108,  1954. 


90 


13.  Rotta,  J.  C.,  "Similar  Solutions  of  Turbulent  Boundary  Layers,"  Jn.  Aero. 
So.,  Vol . 22,  pp.  215-216,  1955. 


14.  Townsend,  A.  A.,  "The  Properties  of  Equilibrium  Boundary  Layers,"  Jn. 
Fluid  Nlech . , Vol.  1,  pp.  561-573,  1956. 

15.  Coles,  D.,  "Remarks  on  the  Equilibrium  Turbulent  Boundary  Layer," 

Jn.  Aero.  Sci.,  Vol.  24,  pp.  495-506,  1957 

16.  Rotta,  J.  C. , "Turbulent  Boundary  Layers  in  Incompressible  Flow," 

Prog,  in  Aero.  Sci.,  Vol.  2,  pp.  1-219,  Pergamon  Press,  1962. 

17.  Perry,  A.  E. , Schofield,  W.  H. , and  Joubert,  P.  H.,  "Rough  Wall  Turbulent 
Boundary  Layers,"  Jn.  Fluid  Mech.,  Vol.  37,  pp.  383-413,  1969. 

18.  Schlichting,  H. , Boundary  Layer  'Iheory , 6th  Edition,  McGraw-Hill  Book 
Co.,  Inc.,  New  York,  1968'. 

19.  Andersen,  P.  S.,  Kays,  W.  M. , and  Moffat,  R.  J.,  "The  Turbulent  Boundary 

Layer  on  a Porous  Plate:  An  Experimental  Study  of  the  Fluid  Mechanics 

for  Adverse  Free-Stream  Pressure  Gradients,"  Report  No.  HMT-15, 
Thermosciences  Division,  Dept,  of  Mech.  Eng.,  Stanford  University,  1972. 

20.  Klebanoff,  P.  S.  "Characteristics  of  Turbulence  in  a Boundary  Layer 
with  Zero  Pressure  Gradient,"  NACA  Report  1247,  1955. 

21.  Blake,  W.  K. , "Turbulent  Boundary- Layer  Wall-Pressure  Fluctuations  on 
Smooth  and  Rough  Walls,"  Jn.  Fluid  Mech.,  Vol.  44,  pp.  637-660,  1970. 

22.  Hinze,  J.  0.,  Turbulence  - An  Introduction  to  its  Mechanism  and  Theory, 
McGraw-Hill  Book  Co.,  1959. 

23.  Corrsin,  S.,  and  Kistler,  A.  L.,  "The  Free  Stream  Boundaries  of  Turbu- 
lent Flows,"  NACA  Tech.  Note  3133,  January  1954. 

24.  Thielbahr,  W.  H. , Kays,  W.  M.,  and  Moffat,  R.  J.,  "The  Turbulent 

Boundary  Layer:  Experimental  Heat  Transfer  with  Blowing,  Suction  and 

Favorable  Pressure  Gradient,"  Report  No.  FMT-5,  Thermosciences  Division, 
Dept,  of  Mech.  Eng.,  Stanford  University,  1969. 

25.  Kearney,  D.  W. , Moffat,  R.  J.,  and  Kays,  W.  M. , "The  Turbulent  Boundary 

Layer:  Experimental  Heat  Transfer  with  Strong  Favorable  Pressure 

Gradients  and  Blowing,"  Report  No.  HMT-12,  Thermosciences  Division, 

Dept,  of  Mech.  Eng.,  Stanford  University,  1970. 

26.  Moretti,  P.  M. , and  Kays,  W.  M. , "Heat  Transfer  to  a Turbulent  Boundary 
Layer  with  Varying  Free-Stream  Velocity  and  Varying  Surface  Temperature 
An  Experimental  Study,"  Int.  J.  Heat  Mass  Transfer,  Vol.  8,  pp.  1187- 
1202,  1965. 


' 


91 


i 

[ 


I 

I 

I 

I 

I 


I 

I 

I 

I 

I 

I 

I 

! 

I 


27.  .Julien,  H.  L.,  Kays,  W.  M. , and  Moffat,  R.  J.,  "Hie  Turbulent  Layer 

on  a Porous  Plate:  bxper imental  Study  of  the  Effects  of  a Favorable 

Pressure  Gradient,"  Report  No.  UNIT- 4,  Thermosciences  Division,  Dept, 
of  Mech.  ling. , Stanford  University,  1 1)69 . 

28.  Loyd,  R.  J.,  Moffat,  R.  .J.,  and  Kays,  W.  M.,  "The  Turbulent  Bounda^ 

Layer  on  a Porous  Plate:  An  Experimental  Study  of  the  Fluid  Dynamics 

with  Strong  Favorable  Pressure  Gradients  and  Blowing,"  Report  No.  IMT-13, 
Thermosciences  Division,  Dept,  of  Mech.  Eng.,  Stanford  University,  1970. 

29.  Kays,  W.  M. , Convective  lleat  and  Mass  Transfer,  McGraw-Hill  Book  Co., 

New  York,  1 9bTT. 

30.  Monin,  A.  S.,  and  Yaglom,  A.  M. , Statistical  Fluid  Mechanics,  Vol.  1, 

The  MIT  Press , 1971. 

31.  Tennekes,  H. , and  Lumley,  J.  L.,  A First  Course  in  Turbulence,  The  MIT 
Press,  1972. 

32.  Townsend,  A.  A.,  The  Structure  of  Turbulent  Shear  Flow,  Cambridge 
University  Press,  1956. 

33.  Bradshaw,  P. , "The  Turbulence  Structure  of  Equilibrium  Boundary  Layers," 
NPL  Aero  Report  1184,  1966. 

34.  Orlando,  A.  F.,  Moffat,  R.  J.,  and  Kays,  W.  M.,  "Turbulent  Transport  of 
Heat  and  Momentum  in  a Boundary  Layer  Subject  to  Deceleration,  Suction 
and  Variable  Wall  Temperature,"  Report  No.  HMT- 17 , Thermosciences 
Division,  Dept,  of  Mech.  Eng.,  Stanford  University,  1974. 

35.  Kim,  H.  T.,  Kline,  S.  J.,  and  Reynolds,  W.  C.,  "The  Production  of 
Turbulence  Near  a Smooth  Wall  in  a Turbulent  Boundary  Layer,"  Jn.  Fluid 
Mech. , Vol.  50,  pp.  133-160,  1971. 

36.  Kline,  S.  J.,  Reynolds,  W.  C.,  Schraub,  F.  A.,  and  Runstadler,  P.  W. , 

"The  Structure  of  Turbulent  Boundary  Layers,"  Jn.  Fluid  Mech.,  Vol.  30, 
pp.  741-773,  1967. 

37.  Offen,  G.  R. , and  Kline,  S.  J.,  "Combined  Dye-Streak  and  Hydrogen- Bubble 
Visual  Observations  of  a Turbulent  Boundary  Layer,"  Jn.  Fluid  Mech., 

Vol.  62,  pp.  223-239,  1974. 

38.  Grass,  A.  J.,  "Structural  Features  of  Turbulent  Flow  over  Smooth  and 
Rough  Boundaries,"  Jn,  Fluid  Mech. , Vol.  50,  pp.  233-256,  1971. 

39.  Reynolds,  W.  C.,  Kays,  W.  M. , and  Kline,  S.  J.,  "Heat  Transfer  in  the 
Turbulent  Incompressible  Boundary  Layer  with  Arbitrary  Wall  Temperature 
and  Heat  Flux,"  Final  Report  - Part  III,  Dept,  of  Mech.  Eng.,  Stanford 
University,  1957. 

40.  Whitten,  D.  G.,  Kays,  W.  M.,  and  Moffat,  R.  J.,  "The  Turbulent  Boundary 
Layer  on  a Porous  Plate:  Experimental  Heat  Transfer  with  Variable 
Suction,  Blowing,  and  Surface  Temperature,"  Report  No.  HMT-3,  Thermo- 
sciences Division,  Dept,  of  Mech.  Eng.,  Stanford  University,  1967. 


92 


41.  Hottcl,  H.  C.,  and  Kalitinsky,  A.,  "Temperature  Measurements  in  High- 
Velocity  Air  Streams,"  Jn.  Appl.  Mech.,  Vol.  12,  pp.  A25-A32,  1945. 

42.  Jorgensen,  F . E. , "Directional  Sensitivity  of  Wire  and  Hot-Film  Probes," 
PISA  Information  No.  11,  1971. 

43.  Repik,  Y.  E.  U.,  and  Ponomareva,  V.  S. , "The  Effect  of  Proximity  of 
Walls  on  the  Readings  of  a Hot-Wire  Anemometer  in  Turbulent  Boundary 
Layers,"  Heat  Transfer  - Soviet  Research,  Vol.  2,  No.  4,  1970. 

44.  l.igrani,  P. , private  communication,  Dept,  of  Mech.  Eng.,  Stanford 
University,  1975. 


l 


l 


APPENDIX  I 

DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 
AND  MEASUREMENT  TECHNIQUES 

The  experimental  apparatus  and  measurement  techniques  used  in  this 
investigation  are  discussed  below.  Since  comprehensive  descriptions  of  the 
design,  construction,  and  basic  qualification  of  the  Roughness  Rig  have  been 
reported  by  Healzer  [1]  and  Pimenta  [2],  these  details  will  not  be  repeated 
here.  The  interested  reader  is  referred  to  the  cited  reports.  Since  some 
differences  exist  between  the  hot  wire  techniques  used  by  Pimenta  [2]  and 
the  present  author,  a more  detailed  discussion  of  this  area  is  presented. 
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I . 1 Experimental  Apparatus 

A schematic  diagram  of  the  Roughness  Rig  was  presented  previously  in 
Figure  1-1  and  a brief  description  was  given  in  Chapter  1.  Details  of  a 
typical  plate  and  casting  assembly  are  shown  in  Figure  1-1  for  reference. 
Each  of  the  24  plates  which  make  up  the  rough  test  surface  has  its  own  plate 
power  control,  transpiration  air  control,  and  five  thermocouples  which  are 
averaged  to  give  an  effective  plate  temperature. 

Stanton  numbers  were  determined  from  an  energy  balance  on  each  plate. 
The  energy  gained  by  the  transpired  air  while  passing  through  the  plate  and 
the  plate  losses  were  subtracted  from  the  measured  plate  power  input  to 
determine  the  wall  heat  flux,  c|".  Stanton  number  was  then  calculated  from 
the  definition 


The  plate  losses  considered  were  the  radiative  losses  from  the  upper 
and  lower  surfaces  of  the  plate,  the  conduction  losses  from  the  plate  to 


the  casting,  and  the  conduction  loss  through  the  stagnant  air  below  the  test 
plate  when  there  was  no  transpiration.  The  same  models  for  these  losses 
were  used  in  this  investigation  as  were  used  by  Healzer  [1]  and  Pimenta  [2]. 
(See  Appendix  II  for  a discussion  of  qualification  tests). 

Based  on  the  results  of  the  qualification  tests,  the  uncertainty  of 
the  Stanton  number  data  is  within  + 0.0001  Stanton  number  units  for  the 
conditions  of  this  investigation. 

Two  modifications  were  made  to  the  Rig  after  completion  of  Pimenta' s [2] 
zero  pressure  gradient  study  and  prior  to  the  present  work.  First,  a 1/2" 
wide,  1/32”  thick  phenolic  trip  was  installed  with  the  front  edge  three 
inches  inside  the  nozzle  exit.  This  was  done  to  insure  stable  conditions 
at  the  beginning  of  the  acceleration  region  from  run  to  run. 

The  second  modification  was  to  the  top  of  the  test  section.  The 
single-section  plexiglass  top  wall  was  replaced  by  a top  wall  of  1/2"  thick 
plexiglass  constructed  in  five  sections  joined  with  plexiglass  inserts. 

This  wall  allowed  more  precise  control  of  the  pressure  gradient  which  was 
set  along  the  tunnel  length.  The  new  top  wall  was  actually  a reworked 
version  of  the  one  used  by  Julien  [28]  and  Thielbahr  [24]  in  their  smooth 
wall  investigations. 

I . 2 Measurement  Techniques 

1.2.1  Pressure  Measurements 

Two  Statham  unbonded  strain  gauge  differential  pressure 
transducers  were  used  - a model  FM-5  with  a 0 to  0.5  psi  range  and  a mod^l 
IW-97  with  a 0 to  0.05  psi  range.  Each  unit  had  a zeroing  circuit  and  was 
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calibrated  at  the  beginning  of  this  investigation  using  a 30"  Meriam 
Micromanometer  (Model  34I-B2)  as  a standard.  The  calibrations  were  linear 
for  10%  to  80%  of  full  scale  for  both  units  and  were  stable  to  within  0.001 
inches  of  water.  Signals  from  the  transducers  were  integrated  for  ten 
seconds  using  a Hewlett-Packard  Model  2401C  integrating  digital  voltmeter 
(IDVM)  with  an  external  quartz  crystal  oscillator  clock  which  provided  1,  10, 
and  100  second  integrating  period  options. 

The  transducers  were  used  to  read  the  tunnel  static  pressures  from  the 

0.040"  diameter  taps  located  2"  apart  in  the  flow  direction  on  the  tunnel 
side  wall.  This  gave  a pressure  reading  at  the  front,  middle,  and  rear  of 
each  4"  plate  length.  The  value  of  the  pressure  gradient  was  calculated 
from  a quadratic  fit  to  the  pressures  measured  at  three  adjacent  taps. 

Total  pressures  in  the  freestream  were  measured  with  a Kiel  probe 
located  in  the  potential  flow  region  using  the  pressure  transducers. 

1.2.2  Temperature  Measurements 

Mean  temperature  profiles  were  measured  with  a 0.003"  diameter, 
butt-welded  Chromel-constantan  thermocouple  probe  mounted  in  a traversing 
probe  holder  similar  to  that  used  with  the  horizontal  hot  wire  (see  Section 
1.2.3).  The  probe  was  designed  with  a length  of  approximately  0.625"  to 
minimize  conduction  errors  (see  Blackwell's  [9]  analysis),  and  was  calibrated 
in  an  oil  bath  using  a Hewlett-Packard  Model  2801A  Quartz  Thermometer  as  a 
standard.  Recovery  factor  for  this  probe  was  assumed  to  be  0.66  based  on 
the  work  of  Hottel  and  Kalitinsky  [41]. 

The  freestream  temperature  (for  use  in  Stanton  number  runs  and  monitoring 
during  hot  wire  data  acquisition)  was  measured  with  a probe  made  of  0.004 
inch  iron-constantan  wire  welded  into  a bead.  This  probe  was  also  calibrated 

against  the  quartz  thermometer  and  a recovery  factor  of  0.86  was  used. 
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Estimated  accuracy  of  temperature  measurements:  + 0.15  °F. 


E 

1.2.3  Hot  Wire  Measurements 

Two  hot  wire  probes  were  used  in  this  investigation:  (1)  a 

DISA  55P05  horizontal,  boundary- layer- type  probe  of  5 micron  tungsten  wire 
with  gold  plated  ends,  and  (2)  a DISA  55F02  45°  slant  probe  of  5 micron 
tungsten  wire  with  gold  plated  ends.  The  wires  were  mounted  in  the  same 
probe  traversing  mechanisms  described  by  Pimenta  [2]  and  which  are  shown  in 
Figures  1-2  and  1-3  for  reference.  The  slant  wire  could  be  rotated  about 
the  probe  axis  with  stops  positioned  45°  apart  (6-  = mr/4,  where  n = 0 7). 

Two  DISA  5SM01  anemometers  with  CTA  Standard  Bridges  were  operated  in 
the  constant -temperature  mode.  Each  of  the  two  hot  wires  was  paired  with 
an  anemometer.  The  anemometer  output  voltages  were  linearized  using  two 
TSI  Model  1072  fourth  order  linearizers. 

A DISA  Model  55D15  true  rms  meter  was  used  to  determine  the  mean 
square  values  of  the  fluctuating  voltages.  The  rms  meter  was  calibrated 
against  standard  sine  waves  with  known  rms  values.  Resulting  accuracy  of 
the  rms  meter  output  was  1%  of  the  measured  value. 

Mean  velocities  were  determined  by  integrating  the  linearizer  output 
for  10  seconds  with  the  IDVM,  while  mean  square  values  were  determined 
by  integrating  the  true  rms  meter  output  for  100  seconds  with  the  IDVM. 

The  rms  meter  was  used  with  a 10  second  time  constant  setting,  and  four  time 
constants  (40  seconds)  were  allowed  to  elapse  before  the  100  second  integra- 
tion was  begun. 

The  hot  wire  probes  were  calibrated  using  the  calibrator  described  by 
Pimenta.  This  consists  of  a length  of  3"  diameter  PVC  pipe  with  flow 
straighteners  and  screens  at  its  inlet  and  is  followed  by  a 20:1  contraction 

• I 
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ASMF;  nozzle.  The  probes  were  placed  in  the  free  jet  at  the  exit  of  the 
nozzle  where  the  velocity  is  uniform  across  the  central  region  of  the  jet. 
The  air  temperature  was  maintained  constant  and  at  the  same  temperature  at 
which  the  data  were  taken  in  the  tunnel  (~  68°F). 

The  directional  sensitivity  of  a hot  wire  and  the  resulting  equations 
have  been  widely  reported  in  the  literature  and  will  be  only  briefly  covered 
here.  According  to  Jorgensen  [42],  the  directional  sensitivity  of  a hot 
wire  can  be  written  as 


2 ,22.22 
u2  + kl  v2  + k2  W2 


(1.2) 


where  V£,  and  are  the  velocity  components  in  the  wire  coordinate 
system  (Figure  1-4)  and  and  k£  are  constants  which  depend  on  wire  and 
prong  construction  characteristics.  For  the  DISA  55F02,  the  values  are 
taken  as  k^  = 0.20  and  k2  = 1.02. 

Equation  (1.2)  can  be  rewritten  in  terms  of  u^,  v^,  w^,  the  velocity 
components  in  the  laboratory  coordinate  system,  as 

ueff  = ^1  + ®vl  + ^W1  + ^uivi  + ^vlwl  + ^ulwl  (1.3) 

where 

2 2 2 
A = cos  <J>  + k^  sin  4> 

B = (sin^<j>  + k^  cos^4>)  cos^9  + k^  sin^O 

C = (sin^4>  + k^  cos^4>)  sin^O  + k^  cos^O 

D = (1  - k^)  sin  2<f)  cos9 

E = (sin^4>  + k^  cos^9  - k^)  sin  29 

F = (1  - k^)  sin  24>  sin9 
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For  a boundary  layer  flow  where 


= U + u' 


v^  = v' 


w^  = w' 


it  can  be  shown  [2,34]  that 


u'L  = A u'2  + 5*-  v’2  + L-  w'2  + D u'v’  + Sr  v’w’  + F u’w'  + 0(3)  (1.4) 


Ueff  = ^ u + °(2> 


(1.5) 


In  a two-dimensional  boundary  layer,  the  v'w'  and  u'w'  terms  can  be 

assumed  zero  by  symmetry.  In  the  present  tests,  it  was  verified  that  both 

these  terms  were  essentially  zero  to  within  the  accuracy  of  the  rms  meter 

(see  Appendix  II).  With  this  information,  U and  u'2  can  be  measured  with 

the  horizontal  wire,  and  measurements  at  three  rotations  of  the  slant  wire 

— 2 — 1 

allow  determination  of  v'  , w'  , and  u’v' . In  this  investigation,  slant 

wire  measurements  were  taken  at  O'  = 45°,  90°,  and  135°.  These  angles  were 

— 2 

chosen  after  an  initial  investigation  to  determine  typical  values  of  e'  at 

— t — j 

different  0-' s.  Solution  for  v'  , w'  , and  u'v'  involves  finding  small 
differences  between  large  numbers.  The  choice  of  0-  = 45°,  90°,  and  135°  was 
an  attempt  to  maximize  this  difference  and  thus  minimize  the  error  involved. 
In  addition,  these  values  of  0- minimize  the  effect  of  the  velocity  gradient 
in  the  0°  - 180°  plane. 

Estimated  uncertainties  in  the  indicated  quantities  are: 
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1 POROUS  PLATE 

2 HEATER  WIRES 

3 THERMOCOUPLES 

4 SUPPORT  WEB 

5 TRANS.  AIR 
THERMOCOUPLE 

6 HONEY  COMB 

7 CAST  WATER 
PASSAGE 

8 PRE- PLATE 

9 TRANS  AIR 
DEFLECTOR 

10  CASTING 

11  CASTING 
THERMOCOUPLE 

12  TRANS  AIR 
DELIVERY  TUBE 


Figure  1-1.  Cross-Section  of  Typical  Test  Plate- -Casting  Configuration 


DETAIL  OF 
WALL  STOP 


Figure  1-2.  Schematic  of  the  Horizontal  Hot-Wire  Probe  Configuration 
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SPINDLE 


LOCK  DRUM 


Figure  1-3.  Schematic  of  the  Rotatable  Slant  Hot-Wire  Probe  Configuration 
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Figure  1-4.  Coordinate  Systems  for  Analysis  of  Slant  Wire 
Directional  Characteristics 
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QUALIFICATION  TESTS 

11 .1  Stanton  Number  Data 

The  original  energy  balance  tests  on  the  Roughness  Rig  were  performed 
by  Healzer  [1],  and  the  St  data  of  Pimenta  [2]  were  taken  shortly  after  these 
tests.  Since  the  present  St  data  for  K = 0,  F = 0 and  0.0039  are  in 
excellent  agreement  with  those  of  Pimenta  (see  Figure  3-1)  for  large  values 
of  enthalpy  thickness,  it  was  concluded  that  the  models  for  the  energy  losses 
were  still  valid. 

This  conclusion  was  confirmed  when  both  blown  and  unblown  energy 
balance  tests  were  conducted  several  months  after  this  investigation  was 
completed.  Results  of  these  tests  [44]  were  in  excellent  agreement  with 
the  original  tests  of  Healzer. 

11 . 2 Two  Dimensionality  Check 

After  installation  of  the  boundary  layer  trip,  a zero  pressure 
gradient  flow  was  established  in  the  test  section  and  velocity  profiles  were 
taken  at  x = 24"  both  on  the  centerline  and  3 inches  on  either  side 
(z  = 0,  l 3 inches).  Across  this  center  6"  section  of  the  tunnel,  the 
variation  of  U was  less  than  1%  and  the  variation  of  momentum  thickness 

oo 

was  less  than  1%. 

, When  a favorable  pressure  gradient  (acceleration)  is  imposed  on  a 

flow  in  a test  section  of  finite  size,  a divergence  of  the  streamlines  is 
expected  due  to  the  thinning  of  the  boundary  layers  on  the  smooth  side  and 
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top  walls.  Since  the  accelerations  applied  in  this  investigation  were 
milder  (lower  values  of  K)  than  those  studied  in  the  smooth  wall  layer 
investigations  of  Julien  [27],  Thielbahr  [24],  Loyd  [28],  and  Kearney  [25], 
the  three-dimensionality  induced  by  acceleration  was  anticipated  to  be 
negligible  in  this  study.  To  check  this  effect  on  the  turbulence  measure- 
ments, the  values  of  v'w'  and  u'w'  were  determined  for  each  acceleration 
condition  at  the  x-position  where  turbulence  profiles  were  taken.  In 
every  case,  v'w'  and  u'w'  were  essentially  zero  within  the  accuracy  of  the 
rms  meter,  and  in  no  case  were  they  greater  than  2%  of  u’ v'".  It  was 
concluded  that  the  accelerated  boundary  layers  of  this  investigation  could 
be  considered  two-dimensional  with  negligible  loss  of  accuracy. 


II . 3 Sensitivity  of  Calculated  Data  to  Origin  of  Wall 

The  sensitivity  to  the  assumed  origin  of  the  wall  of  all  accelerated 


data  calculated  by  integration  of  profiles  outward  from  the  wall  was 


considered.  All  calculated  quantities  reported  in  this  thesis  were  deter- 
mined assuming  y = 0 at  the  crests  of  the  spherical  roughness  elements. 

Additional  calculations  were  made  considering  the  origin  of  the  wall  to 
be  0.006  inches  below  the  crests  of  the  roughness  elements.  Typical  results, 
presented  as  the  percentage  change  due  to  the  0.006"  wall  shift,  were: 

6 - n 

61  + M 

62  - 1% 

H + 4% 

A -*■  1% 
a2  - U 

Pr^  ■*  3% 

Cf/2  - 2%  . 
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II .4  Hot  Wire  Measurements 

Great  care  was  taken  when  hot  wire  measurements  were  being  made  to 
insure  that  the  wire  calibration  and  instrument  calibrations  were  maintained. 

Periodic  checks  of  wire  calibration  were  made  in  the  freestream  by  comparison 
with  measurements  of  obtained  with  a Pitot  probe  and  the  calibrated 
pressure  transducers.  The  calibration  settings  of  the  linearizers,  rms 
meter,  and  IDVM  were  checked  after  every  one  or  two  profiles  to  insure  that 
these  instruments  had  a minimum  of  drift.  The  temperature  of  the  flow  was 
maintained  constant  within  + 0.5°F  during  measurement  periods  by  monitoring 
the  output  of  a thermocouple  in  the  freestream  on  an  auxiliary  digital 
voltmeter. 

] 
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INDUCED  TRANSPIRATION  EFFECTS  ON  ACCELERATION  DATA 

Imposition  of  a pressure  gradient,  dP/dx,  along  the  test  section  results 
in  an  induced  flow  through  the  porous  plates  for  a nominal  F = 0 condition. 

When  there  is  blowing  through  the  plates,  the  nominal  distribution  of  F 
along  the  test  section  is  altered  by  the  flow  induced  by  the  pressure  gradient. 
This  section  presents  an  analysis  which  quantifies  these  effects. 

A cross-section  diagram  of  a typical  plate  and  casting  assembly  was 
shown  in  Figure  1-1.  In  the  following  analysis  it  will  be  assumed  that,  for 
F = 0,  the  preplate  is  impermeable  since  its  porosity  is  much  less  than  that 
of  the  test  plate.  Static  pressures  on  the  front  edge,  center,  and  rear 
edge  of  the  test  plate  will  be  denoted  by  P^,  and  P^,  respectively.  For 
a linear  variation  of  pressure  along  the  plate,  the  pressure  in  the  plenum 
between  the  preplate  and  test  plate  is  assumed  to  be  (P^  + P^)/2  for  F = 0. 
Under  the  favorable  pressure  gradients  of  this  study,  for  a nominal  F = 0 
condition  there  is  suction  on  the  front  half  of  the  plate  and  blowing  on  the 
rear  half. 

The  pressure  drop  versus  flow  rate  characteristics  of  a typical  plate 
were  determined  experimentally.  The  pressure  differential  across  a test 
plate  was  measured  for  11  different  settings  of  the  transpiration  control 
valve  (while  there  was  no  mainstream  flow).  It  was  found  that  the  data 
» follow  the  relationship 


m = 25.58  AP0,942  (III.l) 
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is  the  pressure  drop  across  the  test  plate  (inches  F^O) . This  expression  is 
valid  for  AP  from  0.019  to  0.56"  I^O  and  m from  0.58  to  14.68  ft3/min. 

Calculations  were  made  using  Equation  (1 1 1 . 1) , the  assumptions  stated 
above,  and  assuming  negligible  induced  flow  axially  in  the  test  plate  due 
to  the  longer  flow  path  in  that  direction  than  in  the  direction  normal  to 
the  plate  surface.  Results  are  shown  below,  where  F • and  F are  the 
blowing  fractions  at  the  leading  and  trailing  edges  of  the  plate,  respectively. 
Maximum  induced  suction  occurred  at  the  leading  edge,  maximum  induced  blowing 
at  the  trailing  edge,  and  the  value  of  F induced  at  the  middle  of  the  plate 
was  zero. 


Run 

C*-xa) 

p 

nominal 

F 

max 

Fmin 

K - 0.15  x 10'J 
r 

22 

42 

0.0 

0.0 

0.0005 

0.0006 

-0.0006 

-0.0006 

K = 0.29  x 10'3 
r 

10 

22 

0.0 

0.0 

0.0010 

0.0011 

-0.0010 

-0.0012 

Kr  = 0.29  x 10'3 

10 

22 

0.0039 

0.0039 

0.0046 

0.0049 

-0.0031 

-0.0030 

K = 0.28  x 10'6 

10 

22 

0.0 

0.0 

0.0013 

0.0019 

-0.0014 

-0.0021 
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Tabulation  of  Experimental  Data 

This  appendix  provides  tabular  listings  of  the  experimental  data  of 
this  investigation.  Data  are  presented  in  the  following  order:  (1)  Stanton 

numbers,  (2)  mean  temperature  profiles,  (3)  mean  velocity  profiles,  and 
(4)  Reynolds  stress  tensor  component  profiles. 

Abbreviations  used  in  the  listings  are: 


RUN 

Date  (month,  day,  year) 

TINFO 

Freestream  total  temperature 

(°F) 

TDB 

Dry  bulb  temperature 

(°F) 

TWB 

Wet  bulb  temperature 

(°F) 

PL 

Plate  number 

-- 

X 

Axial  position  from  nozzle  exit 

(In.) 

ST 

Stanton  number 

-- 

DEH2 

Enthalpy  thickness, 

(In.) 

REDEH2 

Enthalpy  thickness  Reynolds  number,  Re. 

n2 

-- 

F 

Blowing  fraction 

-- 

UINF 

Freestream  velocity 

(ft/sec) 

TW 

Wall  temperature 

(°F) 

TINE 

Freestream  static  temperature 

(°F) 

KR 

Fully  rough  acceleration  parameter,  Kj. 

-- 

K 

Smooth  wall  acceleration  parameter,  K 

-- 

DEH 

DE 

DE2 

PT 

Y 

T 

TEAR 

CE/2 

UTAU 

DEI 

H 

G 

BETA 
REDE  2 
REK 
U 

UP2/UI2 

VP2/UI2 

WP2/UI2 

-UV/UI2 

Q2/UI2 

RUV 

RQ2 


Thermal  boundary  layer  thickness,  A 
Velocity  boundary  layer  thickness,  6 
Momentum  thickness,  62 
Point  number  in  profile 
Distance  normal  to  test  surface 
Mean  static  temperature 


(T  - T)/(T  - T ) 

Skin  friction  coefficient,  C^/2 
Friction  velocity,  Ux 
Displacement  thickness,  6^ 

Shape  factor,  6^/62 
Clauser  shape  factor 

Pressure  gradient  parameter,  B = C51/Tw) (dP/dx) 

Momentum  thickness  Reynolds  number,  Re* 

°2 

Roughness  Reynolds  number,  Re^ 


Mean  velocity 
v'2/U2 


-u'v/u: 

q2/u2 

-u'v'/q2 


.2  ,,1 2 


(In.) 

(In.) 

(In.) 

(In.) 

(°F) 


(ft/sec) 

(In.) 


(ft/sec) 
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STANTON  NO.  RUN  - KR*0,  F*0 
RUN  = 040975  TOfl  = 70.00 

TINFC  r 66.60  TWR  = 57.00 


PL 

X 

ST 

OEM2 

REOFM? 

F 

UINF 

TM 

TINF 

1 

2 

.00531 

.011 

48  4. 

0.0000 

98.08 

92.2 

66.0 

0. 

2 

6 

.00372 

.029 

1309  . 

0 .0000 

98.08 

92.3 

66.0 

0. 

3 

10 

.00327 

.043 

1945. 

0.0000 

98.08 

92.4 

66.0 

0. 

4 

14 

.32301 

.055 

2517. 

o.ococ 

98.08 

92.4 

66.0 

0. 

5 

18 

.00285 

.967 

3052  . 

C. 0000 

88.08 

92.3 

66.0 

0. 

6 

22 

.00263 

.079 

3552. 

0.0000 

89.0  9 

92.3 

66.0 

0. 

7 

26 

.00263 

.0  88 

4032  . 

0.0000 

88.08 

92 .3 

66 . 0 

0. 

8 

30 

.00248 

.099 

4498. 

0.0000 

88.08 

92.4 

66.0 

0. 

9 

34 

.00242 

.109 

494  5. 

0.0000 

89.08 

92.3 

66.0 

0. 

10 

38 

.00238 

.118 

5384. 

0.0000 

88.08 

92.4 

66.0 

0. 

11 

42 

.00238 

.129 

5818. 

0.0000 

88.08 

92.3 

66.0 

0. 

12 

46 

.00236 

.137 

6249. 

0.0000 

88.08 

92.2 

66.0 

0. 

13 

50 

.00233 

.147 

6677  . 

0.0000 

88.09 

92.3 

66.0 

0. 

14 

54 

.00225 

.156 

7095. 

0. 0000 

88.08 

9?  .4 

66.0 

0. 

15 

58 

.00222 

.165 

7502  . 

0.0000 

88.0  8 

92.5 

66.0 

0. 

16 

62 

.00221 

.173 

7907. 

0.0000 

88.08 

92.4 

66.0 

0. 

17 

66 

.0022? 

.182 

8311  . 

0.0000 

88.  08 

92.5 

66.0 

3. 

18 

70 

.00219 

.191 

8713. 

0.0000 

98.08 

92.5 

66.0 

0. 

19 

74 

. 00218 

.200 

9112. 

0.0000 

88.08 

92.4 

66.0 

0. 

20 

73 

.00214 

.209 

9506  . 

0.0000 

88.08 

92.5 

66  . O 

0. 

21 

82 

.00210 

.217 

9893. 

0.0000 

88.08 

92.4 

66.0 

0. 

22 

86 

.00212 

.225 

10277. 

0.0000 

88.09 

92.5 

66.0 

0. 

23 

90 

.00212 

.234 

10663. 

0.0003 

88.08 

92.4 

66.0 

3. 

STANTON  NO.  RUN  - KR=0,  F*0  - FIRST  6 PLATES  UNHEATEO 
RUN  * 040975  TO0  = 70.00 

TINFC  = 67.50  TMB  = 57.00 


PL 

X 

ST 

OEM  2 

REOEM2 

F 

UINF 

TM 

TINF 

1 

2 

0.00000 

0.000 

0 . 

3.0000 

88.16 

67.9 

66.9 

0. 

2 

6 

0.00000 

0 .000 

0 . 

0.0000 

88.16 

67.7 

66.9 

0. 

3 

10 

0. 00000 

0.000 

0. 

8.0000 

88.16 

67.5 

66.9 

0. 

4 

14 

0.00000 

0.000 

0. 

0.0000 

88.16 

67.3 

66.9 

0. 

5 

18 

0.00000 

0.000 

0 . 

0.0000 

88.16 

67.5 

66.9 

0. 

6 

2? 

0.00000 

0.000 

0. 

0.0000 

88.16 

68.4 

66.9 

0. 

7 

26 

.00455 

.009 

414. 

0.0000 

88.16 

93.3 

66.9 

0. 

8 

30 

.00335 

.025 

1133. 

0.0000 

88.16 

93.3 

66.9 

0. 

9 

34 

.00314 

.038 

1724. 

0.0000 

88.16 

93.2 

66.9 

0. 

10 

38 

.00297 

.050 

2280  . 

0.0000 

88.16 

93.4 

66.9 

0. 

11 

42 

.00284 

.062 

2808. 

0.0000 

88.16 

93.0 

66.9 

0. 

12 

46 

.00278 

.073 

3319. 

0.0000 

88.16 

93.1 

66.9 

0. 

13 

50 

. 00270 

.084 

3818. 

0.0000 

88.16 

93.0 

66.9 

0. 

14 

54 

.00259 

.095 

4299. 

0.0000 

88.16 

93.1 

66.9 

0. 

15 

58 

.00253 

.105 

4765. 

0.0000 

88.16 

93.2 

66.9 

0. 

16 

62 

. 00250 

.115 

5223. 

0.0000 

88.16 

93.1 

66.9 

0. 

17 

66 

.00247 

.125 

5675. 

0.0000 

88. 16 

93.2 

66.9 

0. 

18 

70 

.00244 

.135 

6121. 

0.8080 

88.16 

93.2 

66.9 

0. 

19 

74 

.00240 

.144 

6561. 

0.0000 

88.16 

93.2 

66.9 

0. 

20 

78 

.00235 

.154 

6993. 

0.0000 

88.16 

93.3 

66.9 

0. 

21 

82 

.00229 

.163 

7416. 

0.0000 

88. 16 

93.4 

66.9 

8. 

22 

86 

.00230 

.172 

7834. 

0.0000 

80.16 

93.4 

66.9 

0. 

23 

90 

.00230 

.181 

8253. 

0.0000 

88.16 

93.4 

66.9 

0. 

I 

I 
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STANTON 

NO.  RUN 

- KR1 0 • 

F=0  - FIRST  12  PLATES  UNMEATEO 

RUN 

2 

040975 

T08 

« 70.00 

TINfO  * 

67.20 

TNB 

= 57.00 

PL 

T 

ST 

OEM  2 

REOEH2 

F 

UINF 

TM 

TINF 

1 

2 

0.00000 

0.000 

0. 

0.0000 

88.13 

67.9 

66.6 

0. 

2 

6 

0.00000 

0.000 

0 . 

0.0000 

88.13 

67.6 

66.6 

0. 

3 

10 

0.00000 

0.000 

0. 

0.  0000 

88.13 

67.3 

66.6 

0. 

4 

14 

0.30000 

0.900 

8. 

0.0000 

88.13 

67.2 

66.6 

0. 

5 

18 

0.00000 

0.009 

0. 

0.  0000 

88.13 

67.4 

66.6 

0. 

6 

22 

3.00000 

O.OOS 

0. 

0.0000 

88.13 

67.3 

66.6 

0. 

7 

26 

0.00000 

0.000 

0 . 

0.0000 

68.13 

67.3 

66.6 

0. 

9 

30 

0.00000 

0.000 

0. 

0.0000 

88.13 

67.2 

66.6 

0. 

9 

34 

0.00000 

0.000 

0 . 

G.0000 

88.13 

67.2 

66.6 

0. 

10 

38 

0.00000 

0.000 

0. 

0.0009 

88.13 

67.1 

66.6 

0. 

11 

42 

0.00000 

0.009 

0. 

0.0000 

88.13 

67.3 

66.6 

0. 

12 

46 

0.00000 

0.000 

0 . 

0.0000 

88.13 

66.1 

66.6 

0. 

13 

50 

.00427 

.009 

389. 

0,0000 

88.13 

93.4 

66.6 

0. 

14 

54 

. 00320 

.024 

1069. 

0.0000 

88.13 

93.4 

66.6 

0. 

15 

58 

.00304 

.036 

16  T 8 . 

0.0000 

88.13 

93.5 

66.6 

0. 

16 

62 

.00292 

.048 

2181. 

0.0000 

88.13 

93.3 

66.6 

0. 

17 

66 

.00283 

.059 

2704. 

0.0003 

88.13 

93.4 

66 .6 

0. 

18 

70 

.03276 

.071 

3212. 

0. 0000 

88.13 

93.3 

66.6 

0. 

19 

74 

.00269 

.081 

379  8 . 

9.0000 

88.13 

33.5 

66.6 

0 . 

20 

7' 

.00260 

.092 

4189  . 

0.0009 

98.  IT 

33.5 

66.6 

0. 

21 

82 

.00252 

.102 

4655. 

0.0000 

86.13 

93.5 

66.6 

0. 

22 

86 

.03251 

.112 

5113. 

0.0003 

88.13 

93.5 

66.6 

0. 

23 

90 

.03250 

.122 

5569. 

0.90C3 

88.13 

33. 4 

66.6 

0 . 

STANTON 

NO.  RUN  - 

KR=3, 

F-C  - LINEAR  TWALL 

. VARIATION 

RUN 

= 

041975 

T09 

= 70.00 

TINFO  = 

64.20 

TMP 

= 56.00 

»L 

X 

ST 

OE  H2 

REOEH2 

F 

UINF 

TW 

TINF 

1 

2 

.00529 

.011 

485  . 

0.0000 

88.05 

132.0 

63.9 

0. 

2 

6 

.00171 

.0  29 

1308  . 

0.0000 

88.05 

102.1 

63.9 

0. 

3 

10 

.00327 

.043 

1947. 

0.0000 

88.05 

102.1 

63.9 

0. 

4 

14 

.00301 

.055 

2538  . 

0.0003 

88.05 

101.8 

63.9 

0. 

5 

18 

.00291 

.068 

3114. 

0.0000 

68.05 

101.2 

63.9 

0. 

6 

22 

.00260 

.080 

3669. 

0.0000 

88.05 

100.5 

63.9 

0. 

7 

26 

.00259 

.092 

4212. 

0.0003 

88.05 

99.9 

63.9 

0. 

8 

30 

.00243 

.104 

4746. 

0.0000 

68.05 

99.2 

63.9 

0. 

9 

34 

.00237 

.116 

5292. 

0.0000 

88.05 

98.5 

63.9 

0. 

10 

38 

.00231 

.127 

5819. 

0.0000 

88.05 

97.9 

63.9 

0. 

11 

42 

.00228 

.139 

6374. 

0. 0000 

88.05 

97.1 

63.9 

0. 

12 

46 

.00225 

.151 

6906. 

0.0000 

88.05 

96.5 

63.9 

0. 

13 

50 

.00223 

.163 

7475. 

0.0000 

88.05 

95.8 

63.9 

0. 

14 

54 

.00211 

.175 

8032. 

o.ocoo 

88.05 

95.1 

63.9 

0. 

15 

58 

.00209 

.188 

8621. 

0.0000 

88.05 

94.4 

63.9 

0. 

16 

62 

.00207 

.201 

9196. 

0.0000 

88.05 

93.7 

63.9 

0. 

17 

66 

.00205 

.214 

9785. 

0.0000 

88.05 

93.1 

63.9 

0. 

18 

70 

.0020? 

.227 

10614. 

0.0000 

88.05 

92.3 

63.9 

0. 

19 

74 

.00  20P 

.261 

11020  . 

0.0000 

88.05 

91.7 

63.9 

0. 

20 

78 

.00194 

.255 

11682  . 

0.0000 

88.05 

91.0 

63.9 

0. 

21 

8? 

.00190 

.270 

12361. 

0.0009 

88.05 

90.3 

63.9 

0. 

22 

• 6 

.00189 

.284 

13013. 

0.0000 

88.05 

89.7 

63.9 

0. 

21 

90 

.00189 

.300 

13727. 

0.0000 

88.05 

89.0 

63.9 

0. 
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STANTON 

NC.  RUN  - 

KR  -Ot 

F = 0 - BILINEAR  ' 

THALL  VARIATION 

RUN 

= 

0410  75 

TOB 

= 70.00 

TINFO  = 

64.  30 

TWR 

* 56.00 

PL 

X 

ST 

06H2 

REOEH2 

F 

UINF 

TW 

T INF 

1 

2 

.00529 

.011 

48  5. 

0.0000 

88.06 

102.2 

64.0 

0. 

2 

6 

.00371 

.029 

1309. 

0.0000 

88.06 

102.2 

64.0 

0. 

3 

10 

.00327 

.043 

1949. 

0.0000 

88.06 

102.2 

64.0 

0. 

4 

lit 

.00303 

.055 

2525. 

0. 0000 

88.06 

102.2 

64.0 

0. 

5 

18 

.00278 

.069 

3146. 

0. 0000 

88.06 

100.9 

64.0 

0. 

6 

22 

.00257 

.083 

3779. 

0.0600 

88.06 

99.4 

64.0 

0. 

7 

26 

.00253 

.396 

4393. 

0.0000 

88.06 

98.0 

64.0 

0. 

8 

30 

.00237 

.110 

5011. 

0. 0000 

88.06 

96.8 

64.0 

a. 

9 

3it 

.00228 

.124 

5677. 

0.0000 

88.06 

95.3 

64.0 

0. 

10 

38 

.00223 

.139 

6343. 

0.0000 

88.06 

94.0 

64.0 

0. 

11 

N? 

.00214 

.155 

7089. 

0.0000 

88.06 

92.5 

64.0 

0. 

12 

46 

.00209 

.171 

7836  . 

0.0000 

88.06 

91.1 

64.0 

0. 

13 

50 

.00205 

. Id  9 

8668  . 

0.0000 

88.06 

89.7 

64.0 

0. 

l<t 

54 

.00187 

.207 

9454. 

0.0000 

88.06 

38.5 

64.0 

0 . 

15 

58 

.00208 

.205 

9 372. 

0.0000 

88.06 

89.7 

64.0 

0. 

16 

62 

.00217 

.202 

9261. 

0.0000 

88.06 

91.1 

64.0 

0. 

17 

66 

.00223 

.2  00 

9161  . 

0. 0000 

88.06 

92.6 

64.0 

0. 

18 

70 

.00228 

.200 

9132. 

0.0000 

88.06 

94.0 

64.0 

0. 

19 

74 

.03228 

.200 

9130. 

0.0000 

88.06 

95.5 

64 .0 

0. 

20 

78 

.03228 

• 2C0 

9166. 

0.0000 

88.06 

96.8 

64.0 

0 . 

21 

8? 

.00227 

.201 

9198. 

0. 0000 

88.06 

98.2 

64.0 

0. 

22 

86 

.CC  231 

.202 

9263  . 

0.0000 

88.06 

99.6 

64.0 

0. 

23 

50 

.00233 

.204 

9 32  8 . 

0.0000 

88.06 

101.0 

64.0 

0. 

STANTON 

NO.  RUN  - 

KR  = 0 , 

F = 0 . 0 0 29 

RUN 

= 

041775 

TOB 

= 71.50 

TINFO  = 

67.  30 

TUB 

= 57.00 

PL 

X 

ST 

OEH2 

REOEH2 

F 

UINF 

T W 

TINF 

1 

2 

.00353 

.015 

677. 

.0039 

87.98 

97.7 

66.9 

0. 

2 

6 

.00213 

.042 

1899  . 

.0039 

67.98 

97.7 

66.9 

0. 

3 

10 

.00162 

.065 

2948. 

.004  0 

87.98 

97.8 

66.9 

0. 

4 

14 

.0014? 

.087 

394  9. 

.0040 

87.98 

97.7 

66.9 

0. 

5 

18 

.00139 

.108 

4918. 

.0040 

87.98 

97.8 

b6 .9 

0. 

6 

22 

.00121 

.130 

5875. 

.0039 

67.98 

97.7 

66.9 

0. 

7 

26 

.00117 

.150 

6796. 

.0039 

87.98 

97.8 

66.9 

0. 

8 

30 

.00105 

.170 

770  1 . 

.0039 

87.98 

97.8 

66.9 

0. 

9 

34 

.00105 

.190 

8614. 

.0040 

87.98 

97  .8 

66.9 

0. 

10 

38 

.00099 

.210 

9534. 

. 0039 

87.98 

97.7 

66.9 

0. 

11 

42 

.00 100 

.230 

10420  . 

.0040 

67.98 

97.7 

66.9 

0. 

12 

46 

.00098 

.2  49 

11293. 

.0039 

87.98 

97.8 

66.9 

0. 

13 

50 

.09095 

.269 

12183. 

.0039 

87.98 

97.8 

66.9 

0. 

14 

54 

.0008? 

.288 

13047. 

.0040 

87.98 

97.6 

66.9 

0. 

15 

58 

.00094 

.307 

13936. 

.00  39 

87.98 

97.8 

66.9 

0. 

16 

6? 

.00094 

. J?6 

14  80  2. 

. 0039 

87.98 

97.8 

66.9 

0. 

17 

66 

.00085 

.345 

15663. 

. 0040 

67.98 

97.9 

66.9 

0. 

18 

70 

.00085 

.366 

16586  . 

.00  39 

67.98 

97.8 

66.9 

0. 

19 

74 

.00078 

.385 

17466. 

.00  39 

87.96 

97.7 

66.9 

0. 

20 

78 

.0098 2 

.405 

18365. 

.0039 

87.98 

97.7 

66.9 

0. 

21 

82 

.00078 

.423 

19182. 

.0040 

87.98 

97.7 

66.9 

0. 

22 

86 

.00081 

.441 

20022. 

.0040 

67.98 

97.0 

66.9 

0. 

23 

90 

. 08081 

.460 

20841 . 

.0039 

87.98 

97.8 

66.9 

0. 
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STANTON 

NO.  PUN 

- KP*0, 

F »0 • 0 1 39 

- F1PST 

6 PLATES 

UNHEATtO 

HUN 

X 

041775 

TOO 

a 71.50 

TINFO  * 

66.90 

rue 

« 57.00 

PL 

X 

ST 

OEH2 

REOEH2 

F 

UINF 

TM 

TXNF 

1 

2 

0.00000 

.008 

356  . 

.0039 

87.95 

68.2 

66.5 

0. 

2 

6 

o.oocco 

.025 

1141. 

.00  39 

87.95 

68.0 

66.5 

0. 

3 

10 

0.00800 

.048 

2156. 

.0040 

87.95 

67.8 

66.5 

0. 

4 

14 

0.00000 

.080 

3648. 

.0040 

87.95 

67.7 

66.5 

0. 

c 

ie 

0.00000 

.063 

2842. 

.0040 

87.95 

68.0 

66.5 

0. 

6 

22 

0.00000 

.049 

2230. 

.0039 

87.95 

68.8 

66.5 

0. 

7 

26 

.00245 

.023 

1049. 

.0039 

87.95 

97.7 

66.5 

0. 

8 

30 

.00160 

.047 

2121. 

. 0039 

87.95 

97.9 

66.5 

0 . 

9 

34 

.00148 

.069 

3114  . 

.0040 

87.95 

97.9 

66.5 

0. 

10 

38 

.00133 

.090 

4085. 

.0039 

87.95 

97.9 

66.5 

0. 

11 

4? 

.00127 

.111 

5055. 

.0040 

87.95 

97.7 

66.5 

0. 

12 

46 

.00121 

.132 

5987. 

. 0039 

87.95 

97.8 

66.5 

0. 

13 

50 

.00116 

.153 

6931  . 

. 00  39 

87.95 

97.7 

66.5 

0. 

14 

54 

.00102 

.173 

7 84  8 . 

.0043 

87.95 

97.7 

66.5 

0. 

15 

58 

.00112 

.193 

8739  . 

. 00  39 

87.95 

97.8 

66.5 

3. 

16 

52 

.00110 

.213 

9663. 

.0039 

87.95 

97.8 

66.5 

0. 

17 

66 

.00101 

. 232 

10538. 

. 00  39 

B7.95 

97.9 

66.5 

0. 

18 

70 

.00096 

.252 

11422. 

. 00  39 

87.95 

97.9 

66.5 

0. 

19 

74 

. 00097 

.272 

12333  . 

. 00  39 

8 7.95 

97.8 

66.5 

0. 

20 

78 

. CC  094 

.290 

13174. 

. 00  39 

87.95 

97.9 

66.5 

0. 

21 

82 

.00089 

.310 

14070. 

.0040 

87.95 

97.9 

66.5 

0. 

22 

86 

.00091 

.329 

14952 . 

. 0040 

87.95 

97.9 

66.5 

0. 

23 

90 

.00091 

.349 

15850  . 

. 0039 

87.95 

97.9 

66.5 

c. 

STANTON 

NO.  PUN  - 

KP  = C , 

F-0.0039 

- BILINEAR  r N ALL 

VARIATION 

RUN 

S 

041775 

tob 

* 71.50 

T X NF  0 = 

67.30 

TMB 

= 57.00 

PL 

X 

ST 

OEM  2 

PE  Ofc  H 2 

F 

UINF 

TN 

TINF 

1 

2 

.00355 

.015 

678  . 

. 0039 

87.98 

37.6 

66.9 

0. 

2 

6 

.00215 

.042 

1905. 

.0039 

87.  98 

97.7 

66.9 

0. 

3 

10 

.00165 

.065 

2957. 

. 0040 

87.98 

97.8 

66.9 

0. 

4 

14 

.00145 

.068 

3967. 

.0040 

97.98 

97.7 

66.9 

0. 

9 

18 

.00136 

.113 

5101. 

. 0040 

87.98 

96.5 

66.9 

0. 

6 

22 

. CO  1 1 6 

.138 

6280  . 

. 0039 

97.98 

95.2 

66.9 

0. 

7 

26 

.00111 

.166 

7544. 

.0039 

9 7.98 

93.8 

66.9 

0. 

8 

30 

.00095 

.195 

8862. 

. 0039 

87.98 

92.4 

66.9 

0. 

9 

34 

.00095 

.226 

10243. 

.0040 

87.98 

91.1 

66.9 

0. 

10 

38 

.00086 

.261 

11846  . 

.0039 

87.98 

99.5 

66.9 

0. 

11 

42 

.00085 

.295 

13366 . 

.0040 

97.98 

68.4 

66.9 

0 . 

12 

46 

.00076 

.337 

15303. 

. 00  39 

87.98 

46.8 

66.9 

0. 

13 

50 

.00082 

.383 

17375. 

.0039 

8T.  98 

45.4 

66.9 

0. 

14 

54 

.00047 

.427 

19352. 

.0040 

87.98 

44.3 

66.9 

0. 

15 

58 

.00086 

.417 

18924. 

.0039 

87.98 

45.5 

66.9 

0. 

16 

62 

.00092 

.408 

18506  . 

.00  39 

97.98 

96 . 8 

66.9 

0. 

17 

66 

.00088 

.401 

18179  . 

.0040 

87.96 

48.1 

66.9 

0. 

18 

70 

.00095 

.396 

17942. 

.0039 

87.99 

49.5 

66.9 

0. 

19 

74 

.00085 

.391 

17733  . 

. 0039 

87.98 

90.9 

66.9 

0. 

20 

78 

.00099 

.389 

176»8 . 

.0039 

87.98 

92.2 

66.9 

0. 

21 

82 

.00093 

.386 

17512. 

.0040 

87.98 

93.8 

66.9 

0. 

22 

86 

.00086 

.388 

17619. 

.0040 

87.98 

95.0 

66.9 

0. 

23 

90 

.00084 

.388 

17607. 

. 0039 

87.98 

96.4 

66.9 

0. 

STANTON 

NO.  RUN  - 

KR  = 0 . 

15E-3,  F= 

0 

RUN 

r 

040176 

TO" 

= 73. CO 

• 

TINFO  = 

67.20 

TN  d 

= 56.00 

PL 

X 

ST 

DEH? 

RE  OF  H? 

F 

UINF 

TN 

TINF 

<R 

1 

2 

.00527 

.011 

482  . 

0.0000 

88.46 

93.6 

66.7 

0. 

2 

6 

.00372 

.029 

1302  . 

0.  uooo 

88.10 

93.6 

66.7 

- . 1 23E-04 

3 

10 

.00329 

.043 

1940. 

0.0000 

88.00 

93 .5 

66.7 

- . 600E-05 

4 

14 

.00304 

.056 

2516. 

0.0000 

88.00 

93.5 

66.7 

• 179E-Q5 

5 

10 

.00287 

.067 

3055. 

0.0000 

88.03 

93.6 

66.7 

• S38E-0  5 

6 

22 

,00268 

.078 

3560  . 

0.0000 

88.12 

93.6 

66.7 

. 335E-05 

7 

26 

.00267 

.089 

4047. 

0.0000 

88.17 

93.5 

66.7 

- • 140E-06 

a 

30 

.00262 

.099 

4520. 

0.0000 

88.11 

93.6 

66.7 

- • 112E-0  5 

5 

34 

.00247 

.109 

4975  . 

C.0000 

80.09 

93.5 

66.7 

• 532E-0  5 

10 

38 

.00244 

.119 

5422. 

0.0000 

88.31 

93.5 

66.7 

. 332E-04 

11 

42 

.00247 

.128 

5872. 

0.0000 

89.13 

93.5 

66.7 

. 10  IE -0  3 

12 

46 

.00250 

.135 

6335. 

0.000  0 

91.05 

93.5 

66.7 

.144E-03 

13 

50 

.00261 

.141 

6813. 

0.0000 

93.32 

93.5 

66.7 

.151E-03 

14 

54 

.00245 

.148 

7297. 

0.0000 

95.53 

93.5 

66.6 

. 143E-03 

15 

58 

.00242 

.154 

7783. 

0. 0000 

97.76 

93.5 

66.6 

.147E-03 

16 

62 

.00242 

.160 

8277  . 

0.0000 

100.02 

93.4 

66.6 

. 142E-03 

17 

66 

.00242 

.166 

878  3. 

0.0000 

102.47 

93.4 

66.5 

.15  3E-03 

IS 

70 

.03243 

.17? 

930  3. 

0.000  0 

10**. 87 

93.5 

66.5 

.143E-03 

19 

74 

. 00242 

.178 

9833. 

0.0000 

107. 31 

93.6 

66.5 

. 141E-03 

20 

70 

. 00  23( 

.183 

10366 . 

0.0000 

109.98 

93.6 

66.4 

•156E-03 

21 

82 

.00232 

.186 

10905  . 

0.0000 

112.68 

93.6 

66.4 

• 14  7F.-0  3 

22 

06 

. CO  232 

.19? 

11450  . 

0.0000 

115.18 

93.6 

66.4 

• 139E-0  3 

23 

90 

. 00237 

.198 

1201 3. 

0. 0000 

117.86 

93 .5 

6b  • 3 

. 165E-0  3 

STANTON  NO.  RUN  - KR=0.15E-3,  F*0  - FIRST  6 PLATFS  UNHEAIEO 
RUN  = 032575  TOO  = 71.00 

T I Nf 0 * 65.10  TUB  = 50.00 


PL 

X 

ST 

Of  H? 

REOEH2 

F 

UINF 

TN 

TINF 

KR 

1 

2 

0.00000 

0.000 

0 . 

0.0000 

88.28 

66 . 7 

65.7 

-.241E-04 

2 

6 

0.08000 

0.000 

0 . 

0.0000 

87.94 

66.4 

65.8 

-.143E-04 

3 

10 

0.00000 

0.000 

0. 

0.0000 

87.79 

66.3 

65.8 

-.583E-05 

4 

14 

0.00000 

0.000 

0 . 

0.0000 

87.75 

66.1 

65.8 

•163E-05 

5 

18 

0.00000 

0.000 

0 . 

0.0000 

87.84 

66.2 

65.8 

• '.'91E-05 

6 

2? 

0.00000 

o.eoo 

0. 

0.0000 

87.90 

67.3 

65.8 

•246E-05 

7 

26 

•C0461 

.009 

417. 

0.0000 

68.00 

92.7 

65.7 

- .27  7E-0  5 

8 

30 

. 00  340 

.025 

1140. 

0.0000 

87.82 

92 .7 

65.8 

-• 33  0E-0  5 

9 

34 

.00319 

.039 

1735. 

0.0000 

87.88 

92.6 

65.6 

•935E-05 

10 

38 

. 00303 

.051 

2297. 

0.0000 

88.08 

12. 8 

65.7 

. 32 0E- 0 ■» 

11 

42 

.00297 

.062 

2843. 

0.0000 

88.87 

92.7 

65.7 

.103E-03 

12 

46 

.00296 

.073 

3389. 

0.0000 

90.88 

92.7 

65.7 

. 142E-03 

13 

50 

.00292 

.083 

3944. 

0.0000 

92.99 

92.7 

65.7 

• 149E-0  3 

14 

54 

.00282 

.092 

4499. 

o.oooe 

95.28 

92.8 

65.7 

.145E-03 

15 

58 

.00275 

.101 

5050  . 

0.0000 

97.45 

92.8 

65.6 

.147E-03 

16 

62 

.00273 

.110 

5605. 

0.0000 

99.  73 

92.6 

65.6 

•143E-03 

17 

66 

.00269 

.118 

6167. 

0.0000 

102.13 

92.7 

65.6 

. 152E-0  3 

18 

70 

.00268 

.126 

6737. 

0.0000 

104.62 

92.7 

65.5 

.146E-03 

19 

74 

.00264 

.133 

7315  . 

0.0000 

107.08 

92.7 

65.5 

• 145E-0  3 

20 

78 

. C0258 

.140 

7895. 

0.0000 

109.77 

92.7 

65.4 

• 155E-0  3 

21 

82 

.09253 

.147 

8476. 

0.0000 

112.37 

92.8 

65.4 

. 143E-03 

22 

86 

.00252 

.154 

9063  . 

0.8000 

114.91 

92.7 

65.4 

.137E-03 

23 

90 

.00254 

.161 

9665  . 

0.0000 

117.55 

92.8 

65.3 

.161E-03 

116 


STANTCN 

NO.  RUN  - 

KR=0  . 

297 - 3 * F*0 

RUN 

M 

022575 

TOB 

* 73.00 

TINFO  * 

66.10 

TNB 

= 60.00 

PL 

X 

ST 

DFM2 

PEDEM? 

f 

UINF 

TH 

TINF 

KR 

1 

2 

. 00523 

.011 

472. 

0.0000 

67.24 

93.4 

65.7 

C. 

2 

6 

.00370 

.026 

1277. 

0.0000 

86.99 

93.4 

65.7 

-.549E-05 

3 

10 

.00327 

.042 

1906. 

0.0000 

87.08 

93.3 

65.7 

. 312E-05 

4 

14 

.00302 

.055 

2473. 

C.0000 

67.06 

93.3 

65. 7 

-.631E-06 

«> 

IS 

.00287 

.067 

3005. 

0.0000 

87.35 

93.5 

65.7 

. 722E-0  4 

6 

22 

.00278 

.076 

352  3. 

0.0000 

69.51 

93.4 

65.7 

• 238E-0  3 

7 

26 

.00284 

.084 

4055. 

0.0000 

93.53 

93.4 

65.6 

• 2S2E-0  3 

6 

30 

. 00272 

.091 

4606. 

0.0003 

97.94 

93.4 

65.6 

.2 85E-03 

9 

34 

.00271 

.098 

5168  . 

o.ocoo 

102.30 

93.2 

65.5 

• 279E-03 

10 

3S 

.00263 

.104 

5746. 

o.onoo 

106.93 

93.2 

65.4 

•28CE-03 

11 

42 

.01265 

.110 

6343. 

0.0000 

111.97 

93.4 

65.3 

.280E-03 

12 

46 

.00261 

.115 

6966  . 

0. 0000 

117.32 

93.3 

65.3 

.289E-03 

13 

50 

.00260 

.120 

7612. 

0.0000 

122.79 

93.4 

65.2 

.301E-03 

14 

54 

.00244 

.125 

8264. 

0.0000 

128.44 

93.3 

65.1 

• 166E-0  3 

15 

5S 

.00233 

.134 

8897. 

0.0000 

129.19 

93.2 

65.0 

- • 184E - 0 5 

16 

62 

.00231 

.143 

9515. 

0.0000 

128.80 

93.2 

65.1 

-.170E-04 

17 

66 

.00225 

.152 

10121. 

0.0000 

129.02 

93.3 

65.1 

.210E-04 

IS 

70 

.00224 

.161 

10721. 

0.0000 

129.34 

93.3 

65.0 

•179E-04 

19 

74 

.00221 

.169 

11315. 

0.0003 

129.52 

93.4 

65.0 

-.107E-04 

20 

7S 

.00215 

.178 

11898. 

0.0000 

129.52 

93.4 

65.0 

.415E-05 

21 

8? 

.00212 

.187 

12466 . 

0.0003 

129.42 

93.2 

65.0 

-.9246-05 

22 

86 

.03210 

.196 

13031  . 

0.0000 

129.22 

93.3 

65.0 

-•176E-04 

23 

90 

.03214 

.205 

13594. 

0.0000 

120.49 

93.3 

65.1 

-.358E-04 

ST8NT0N 

NO.  RUN  - 

KR  = 0. 

29E-3,  F*  0 

.0039 

RUN 

s 

030875 

TOB 

* 71.00 

TINFO  = 

66.90 

TNB 

= 59.08 

PL 

X 

ST 

OEH2 

REOEH2 

F 

UINF 

TM 

TINF 

KR 

1 

2 

.00356 

.015 

667. 

.0039 

87. OC 

98.0 

66.5 

0. 

2 

6 

.09212 

.04? 

1871  . 

.0039 

86.87 

96.1 

66.5 

-.2446-05 

3 

10 

. 00167 

.065 

2906. 

. 0039 

66.80 

98.2 

66.5 

- • 240E-05 

4 

14 

.00145 

.087 

3 882. 

.0039 

86.60 

98.2 

66.5 

- . 124E-0  4 

5 

18 

.00145 

.108 

484  0. 

.0039 

86.62 

98.0 

66.5 

• 55  2E-0  4 

6 

2? 

.00134 

.126 

5800  . 

.0039 

86.82 

98.1 

66.5 

• 239E-0  3 

7 

26 

.00133 

.141 

6777. 

.0039 

92. 98 

98.2 

66.4 

• 29  6E-0  3 

8 

30 

.00125 

.155 

7789. 

.0039 

97.49 

98.1 

66.3 

. 294E-3  3 

9 

34 

.00129 

.168 

8846. 

.0039 

102.09 

98.1 

66.3 

.284E-03 

10 

30 

.09125 

.181 

995  8. 

.0039 

106.83 

96.1 

66.2 

. 279E-0  3 

11 

42 

.00132 

.193 

11125. 

.0039 

111.71 

98.0 

66.1 

. 300E-03 

12 

46 

.00125 

.204 

12347. 

. 0039 

117.42 

98.0 

66.0 

. 295E-0  3 

13 

50 

.00122 

.215 

13613. 

.0039 

122.96 

97.9 

65.9 

. 299E-0  3 

14 

54 

.00116 

.225 

14923. 

.0039 

128.7? 

98. 0 

65.8 

, 195E-0  3 

15 

58 

.00110 

.243 

16252  . 

.0039 

130.10 

96.0 

65.8 

. 669E-0  5 

16 

62 

.00103 

.263 

17571. 

. 00  39 

129.95 

90.1 

65.8 

-.158E-04 

17 

66 

.00096 

.282 

18876. 

.0039 

129.61 

96.0 

65.8 

. S54E-05 

IS 

70 

.00086 

.301 

20156. 

.00  39 

129.92 

98.0 

65.8 

- • 197E-0  5 

19 

74 

. 0 0 098 

.321 

21442. 

.0039 

129.84 

96.1 

65.8 

-.  863E-0  5 

20 

78 

.00086 

.340 

22731. 

. 0039 

129. 81 

96.1 

65.8 

- • 29  IE-05 

21 

82 

. 00080 

.359 

23994. 

.0039 

129. 69 

96.1 

65.8 

.62  8E-05 

22 

86 

.00071 

.378 

25236. 

.0039 

129.74 

98.2 

65.8 

- . 526E-0  5 

23 

90 

.00079 

.396 

26476. 

.0039 

129.76 

96.2 

65.8 

.789E-05 

iU 


117 


PL 

X 

ST 

DFH? 

9F0EH2 

F 

UISF 

TW 

T INF 

K 

1 

2 

. 00523 

.011 

470  . 

0.0000 

87.09 

95.1 

66.9 

-.1216-07 

2 

6 

.00369 

.028 

1270  . 

0.0000 

86.90 

95 .1 

66.9 

-.5496-08 

3 

10 

.00326 

.0-2 

1890  . 

0.0000 

86.92 

95.2 

66.9 

-.3696-06 

it 

14 

.00301 

.0  65 

2451  . 

o.ouoo 

86.82 

95.2 

66.9 

-.15*6-07 

c 

1 R 

.00287 

.067 

2 9b  7 . 

0.0000 

86.56 

95.3 

66.9 

.5406-07 

6 

22 

.00279 

.076 

J432. 

0.0000 

88. 7o 

95 .3 

66.8 

.2166-06 

7 

26 

.30288 

.084 

4015. 

0. coco 

93.36 

95.3 

06 .8 

.2726-06 

8 

30 

. CO  27  T 

.090 

4590  . 

0. 0000 

98.66 

95.1 

66.7 

.2786-06 

1 

34 

.30277 

.096 

5171  . 

o.oooc 

1 04.41 

95. 1 

cfc.6 

.2726-06 

10 

3R 

.00270 

.131 

5770. 

0.000  3 

110.80 

95.1 

66.5 

.2766-06 

11 

42 

.00274 

.106 

6416. 

0.30C3 

118.43 

95.2 

66.* 

.2946-06 

1? 

46 

.00270 

.108 

7109. 

0.3000 

127.87 

95.2 

6b. 2 

.2946-06 

13 

50 

.00  27c 

.111 

7 86  3 . 

0.0000 

138.33 

95.2 

66.0 

.2326-06 

14 

54 

.00247 

.113 

868  1 . 

0.0000 

149.90 

95.0 

t5.8 

. 1636-06 

IS 

5R 

.00233 

.121 

9391  . 

0.0000 

151.26 

95.1 

65.8 

-.8456-06 

16 

62 

.00230 

.132 

10116  . 

0.0000 

149.96 

95.1 

65.8 

-.26*6-07 

17 

66 

.00226 

.141 

10776 . 

0.0000 

149. 21 

95.2 

65.8 

-.5496-08 

lfl 

70 

.00225 

.150 

11465  . 

G. OCCO 

149.16 

95.2 

65.8 

.5116-09 

19 

74 

.00221 

.169 

12134. 

0.0003 

149.65 

95.3 

65.8 

.2336-08 

20 

7a 

.00218 

.167 

12792 . 

0.0000 

149.9’ 

95.3 

65.8 

.7386-08 

21 

R? 

.00214 

.175 

13488. 

0.0001 

150.25 

95.2 

65.8 

.1366-07 

22 

86 

.00215 

.l'U 

14167. 

0.0003 

1 c 3 . 8 5 

95.2 

65.8 

. 1336-07 

23 

90 

.00219 

.192 

14839. 

0. 0003 

161.39 

96 .2 

65.6 

.5756-08 

ST8NT0N 

NC.  RUN 

- K=0 . 2R6-6  . F = 0 

- FIRST 

6 PL4T6S 

UNH6  4 T 1 0 

RUN 

- 

0429  75 

TOP 

= 74. CO 

TINFO  ■= 

67.40 

T H G 

= 58.00 

PL 

X 

ST 

06  H? 

R60FH2 

F 

UINF 

TW 

I INF 

K 

1 

2 

o.oooor 

0 .000 

0 . 

0.0000 

87.  1C 

66 .6 

67.0 

-.1216-07 

2 

6 

0.00000 

0 .000 

0 . 

0. 0000 

86.91 

68.4 

67.0 

- . 54  96-0  d 

3 

10 

0.00090 

0.000 

0 . 

0.0000 

86.93 

68.1 

67.0 

-.3896-08 

4 

14 

0.00000 

0.000 

0 . 

0.0000 

86.84 

68.0 

67.0 

-.1546-37 

5 

16 

o.oooco 

0.000 

0 . 

o.oooc 

86.57 

68.2 

67.0 

.5406-07 

6 

22 

n.ooooo 

0 .000 

0. 

o. oooo 

88.78 

69.1 

67.0 

.2166-06 

7 

26 

.00490 

. 0 1G 

471  . 

0.0000 

93.39 

95.0 

67.0 

.2726-06 

8 

30 

.00  169 

.026 

1315. 

0.0000 

98.67 

95.1 

6b. 9 

.2786-06 

9 

34 

.00349 

.0  39 

2075  . 

o.ooog 

104.43 

94.9 

66.8 

.2726-06 

10 

36 

.00326 

.0*0 

2820  . 

0.0000 

110.82 

95.0 

66.7 

.2766-06 

11 

42 

.00321 

.059 

3579. 

0.0000 

118.45 

95.0 

66.6 

.2946-06 

12 

46 

.00313 

.067 
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Reynolds  stress  tensop  components 


KR-0,  F=0  RUN 


RUN 

= 0 1 1 1 TO 

CF/2 

= .00234 

Otl 

* .266 

BETA 

= 0 

.00 

PLATE 

= 17 

UTAU 

= 4.24 

OE2 

* .182 

RE  OE  2 

* 8081. 

X 

= 66 

F 

= 0.3000 

H 

= 1.46 

REK 

s 

67. 

UINF 

= 07.73 

or 

= 1.35? 

G 

= 6.53 

KR 

= 0. 

Y 

Y/OE 

U 

U32/UX2  VP2/UI2 

MP2/UI2 

-UV/UI2 

Q2/UI2 

RUX 

RU2 

.009 

. 007 

2 7.59 

. 0C783 

.020 

.016 

34.  11 

. 00812 

.030 

.022 

37.  63 

. 00041 

.045 

.033 

41.  08 

.00091 

.060 

.0  44 

43.67 

.00916 

.000 

.009 

46.  74 

. 00935 

.100 

.0  74 

40.  06 

.00939 

.130 

.096 

51.32 

.00935 

003C9 

.00505 

.00221 

.01749 

.411 

.126 

.165 

.122 

53.  79 

.00910 

00330 

.00469 

.00227 

.01 709 

.414 

.133 

. ? 0 0 

.140 

56.  28 

.0C864 

00318 

.00466 

.00226 

.01667 

.427 

.136 

.300 

.222 

61.33 

.03601 

00302 

. 00435 

.00217 

.01530 

.441 

.141 

. 400 

.296 

64.86 

.00715 

00305 

.00408 

.00198 

.01428 

.424 

.139 

.600 

.444 

71.48 

. 00572 

00254 

.00337 

.00167 

.01162 

.438 

.144 

.000 

.59? 

76.  92 

. 00433 

00181 

.00256 

.00124 

. 0 0d70 

.443 

.143 

1.100 

.014 

0 3.  74 

.00216 

00105 

. 00113 

.00057 

.CO  434 

.379 

.131 

1.400 

1.036 

07.  21 

.00041 

00036 

.00026 

.00011 

.00103 

.287 

.10  6 

1.000 

1 . 331 

67.  73 

. 00005 

REYNOLDS  STRESS  TENSOR 

COMPONENTS  - KR= 

0 .15E-3, 

F=0.0000 

RUN 

RUN 

= 040375 

CF/2 

= .00266  OE1 

= . 183 

BETA 

= 

.40 

PLATE 

= 17 

UTAU 

- 5 

.19  OE  2 

= .131 

REDE2 

= 6854. 

X 

= 56 

F 

= c.oaoo  h 

= 1.43 

REK 

= 

84. 

UINF 

= 190.71 

OE 

= l. 

119  G 

= 5.58 

KR 

= • 146E 

-03 

r 

Y/OE 

U 

UP2/UI2 

VP2/UI2 

NP2/UI2 

-UV/UI2 

Q2/UI2 

RUW 

RQ2 

.006 

.005 

33.47 

.00737 

.010 

. 009 

37.  68 

.00307 

.020 

.010 

44.09 

. 00873 

.035 

. 031 

50.04 

. 00936 

.045 

.040 

53.06 

.00959 

.060 

.0  54 

5 6.44 

. 09961 

.080 

.071 

5 9.97 

.09967 

.105 

.094 

63.58 

. 00942 

.130 

.116 

66.46 

.00909 

. 00224 

.00366 

.00214 

.01499 

.474 

.143 

.155 

.139 

6 0.  80 

.00073 

.00234 

.00351 

.00210 

.0  1456 

• 4b5 

.144 

.190 

. 170 

71.  75 

. 00818 

.00196 

.00332 

.00196 

.01348 

.487 

.145 

.225 

.201 

74.46 

.00754 

.00202 

. E 0 31 2 

.00186 

.01268 

.477 

.14  7 

.275 

.246 

7 7.47 

. 00679 

.00187 

.00282 

.00173 

.01148 

.486 

.151 

.350 

.313 

80.92 

. 00566 

.00185 

.00263 

.00191 

.01014 

.467 

.149 

.425 

• 380 

04.  16 

. 00401 

.00163 

.00229 

.00131 

.00873 

.468 

. 150 

.500 

.447 

86.79 

. 00410 

.00153 

.00204 

.00115 

.00767 

.459 

.150 

.700 

.626 

92.56 

.00259 

.00108 

.00133 

.00074 

.00500 

.442 

.148 

.900 

.804 

96.95 

. 00143 

.00065 

.00064 

.00038 

.90272 

• 394 

.140 

1. 150 

1.020 

100.06 

. 00032 

! 

148 

I 


t 


I 


1 


REYNOLDS 

STRESS  TENSOR 

COMPONENTS  - KR  = 

0.15E-3, 

FzO.OOOO  RUN 

RUN 

S 

01.0375 

CF/2 

= .00263  OE 1 

* .181 

BETA  * 

.38 

PlATE 

= 

22 

UTAU 

- 5 

.85  OE  2 

= .131 

REDE2  = 7643. 

X 

S 

86 

F 

= 0.0000  H 

* 1.33 

REK  * 

93. 

UlNf 

= 

111.. 07 

OE 

= 1. 

152  G 

= 5.41 

KR  = 

= . 139E 

-0  3 

Y 

V/OE 

U 

UP2/UI2 

VP2/UI2 

MP2/UI2 

-UV/UI2 

Q2/UI2 

RUV 

RQ2 

. 006 

. 005 

40.05 

. 00740 

.008 

. 007 

42.  11 

. 00773 

.018 

.015 

49.  34 

. 00874 

.033 

.028 

5 6.  11 

.03933 

.0<.5 

. 038 

5 9.  86 

. 00957 

.060 

.050 

63.  70 

. OC  962 

.080 

.067 

67.  39 

. OC  963 

.10? 

. 088 

71.58 

. 00939 

.133 

.109 

7 4.  94 

.00907 

.00195 

.00365 

.00202 

.01467 

.480 

.138 

.15? 

. 138 

7 7.49 

. 00869 

.00200 

.03329 

.00191 

.01 39d 

.458 

.137 

.190 

.159 

80.85 

.00803 

.00204 

.00317 

.00194 

.01329 

.455 

.139 

.225 

. 189 

83.61 

. 00750 

.00207 

.00209 

.00176 

.01246 

.44  7 

.141 

.275 

.231 

87.26 

. 00672 

.09195 

.00266 

.00162 

.01133 

. 448 

.143 

.350 

• 294 

91.03 

. 00563 

. 00164 

.00242 

.00138 

.00969 

.954 

.142 

.<*25 

.357 

94.  71 

.00475 

.00155 

.00215 

.00123 

.00845 

.453 

.146 

.500 

.419 

97.  34 

.00397 

.00140 

.00183 

.00105 

.00728 

.431 

.144 

.700 

.587 

103.54 

.00241 

.00099 

.00124 

.00067 

.0  0 464 

.434 

.144 

.900 

.755 

108.02 

.00139 

.00059 

.00065 

.00038 

.00263 

.420 

.144 

1 . 150 

.965 

111.49 

. 00394 

REYNOLDS 

STRESS  TENSOR 

COMPONENTS  - KR  = 0 

. 29E-3  . 

F=0.0000 

RUN 

RUN  = 

322275 

CF/2 

= .03307  OE1 

= .111 

9ETA 

- . 

.41 

PLATE  = 

11 

UTAU 

= 6 

.14  OE  2 

= .077 

REOt  2 

= 4 4 66  . 

X = 

42 

F 

= 0.3000  H 

= 1.43 

REK 

s 

99. 

UINF  = 

110. 7T 

OF 

~ • 

658  G 

= 5.46 

KR 

= . 28  2E 

-03 

Y 

Y/OF 

U 

UP2/UI2 

VP2/UI2 

MP2/UI2 

-UW/UI2 

Q2/UI2 

RUV 

RQ2 

.008 

.012 

43.  70 

. 00770 

.013 

.020 

48.22 

. 00842 

.019 

.029 

52.31 

.00899 

.023 

.035 

54.53 

.00924 

.027 

.041 

56.51 

. 00939 

.031 

. 047 

5 8.05 

. 00953 

.037 

.056 

63.70 

.00958 

.045 

. 068 

63.39 

. 00968 

.055 

.084 

66.28 

.00966 

.065 

.099 

6 8.61 

. 00951 

. 078 

.114 

70.  88 

. OC  939 

.100 

.152 

74.84 

.00885 

. 130 

. 198 

79.48 

.00816 

.00241 

.0038? 

.00209 

.01 438 

.472 

.145 

. 150 

.228 

82.01 

. 0C767 

. 00229 

.00357 

.00192 

.01353 

• 458 

. 142 

.175 

.266 

85.  00 

. 00 ’12 

.00210 

.00318 

.00183 

.01248 

• 464 

.147 

.200 

.304 

87.32 

. 0C648 

.00204 

.00304 

.OOlofl 

.01156 

.462 

.145 

.225 

.342 

89.  52 

.00593 

.00182 

.00273 

.00151 

.01048 

.460 

.14o 

.275 

.418 

93.38 

. OP 4 *7 

.00154 

.00230 

.00129 

.00871 

.471 

.148 

.350 

.532 

98.32 

.00359 

.00117 

.00174 

.00046 

.00650 

.468 

.i«»e 

.450 

.684 

103.65 

.00224 

.00079 

.00109 

. 00061 

.00412 

.459 

.148 

.600 

.912 

108.64 

.00C7’ 

.00039 

.00041 

.00020 

.00157 

.365 

.127 

.750  1 

.140 

110.60 

. 00C 15 

149 


1 


r 


REYNOLDS  STRESS  TENSOR 

C CMPONE  NTS  - KR*0 

• 29c - 3* 

F = Q . 00  39 

RUN 

RUN 

= 030675 

CF/2 

= .00161  OE1 

* .197 

SETA 

s 

-1 

.44 

PLATE 

= 11 

UTAU 

= 4 

.44  OE2 

= .131 

REDE2 

= 

7524. 

X 

= 4? 

F 

= .0039  H 

= 1.50 

RFK 

s 

72. 

UINF 

= 110.67 

DE 

- 1. 

019  r. 

* 8.35 

KR 

= 

• 294E 

-03 

Y 

Y/OE 

U 

UP2/UI2 

VP2/UI2 

MP2/UI2 

-UV/UI2 

Q2/UI2 

RUV 

RQ2 

.008 

. 008 

33.06 

. 00836 

.013 

.013 

36.92 

.00889 

.023 

.023 

42. 63 

. 00969 

.0  32 

. 031 

4 6.28 

. 01031 

.046 

.045 

5 0.  88 

.01088 

.065 

. 064 

56.16 

.01135 

.075 

.074 

58.21 

.01141 

.087 

.085 

60.  64 

.01147 

.100 

.098 

62.89 

.01145 

.115 

. 113 

65.44 

.01136 

.130 

.128 

67.73 

. 01117 

.00318 

.00525 

.00271 

.01960 

.455 

.138 

.150 

. 147 

70.  15 

.01091 

.00301 

.00513 

.00269 

.01905 

.469 

.141 

.195 

. 91 

74. 68 

.01000 

.00304 

.00463 

.00261 

.01767 

.473 

.148 

. 295 

.289 

82.75 

. 00813 

.00270 

.00381 

.00224 

.01464 

.478 

.153 

. 390 

. 383 

88.60 

. 00667 

.00215 

•C0314 

.00186 

.01196 

.491 

.15  6 

.550 

.540 

9 6.44 

.00452 

.00167 

.00233 

.00133 

.00852 

.484 

.156 

. 700 

.687 

102. 38 

.00297 

.00118 

.00158 

.00089 

.00573 

.475 

.155 

.900 

. 883 

107. 80 

.00122 

.00070 

. 00C65 

.00037 

.00257 

.400 

.144 

1.  IOC 

1.0  79 

110. 24 

. 0 0 C 2* 

.03024 

.30015 

. 000C7 

.00063 

.292 

.111 

REYNOLDS  STRESS  TENSOR 

COMPONENTS  - K=0 

.28E-6,  F 

=0.0000 

RUN 

RUN 

= 050375 

CF/2 

= .00305  DEI 

s .098 

BETA 

= 

- 

.57 

PLATE 

= 11 

UTAU 

= 6.42  OE  2 

= .C69 

RE0E2 

X 

4189. 

X 

* 42 

F 

= C.0000  H 

* 1.42 

REK 

X 

104. 

UINF 

= 116.23 

DE 

- .610  G 

* 5.33 

K 

3 

. 29  2E 

-06 

V 

Y/OE 

U 

UP2/UI2  VP2/UI2 

WP2/UI2 

-UV/UI2 

Q2/UI2 

RUV 

RQ2 

.006 

.010 

45.74 

.00763 

.010 

.016 

50.  22 

. 00801 

.017 

.028 

5 6.  04 

.00365 

.028 

.046 

52.70 

. 00922 

.035 

.057 

65.95 

. 00932 

.045 

. 074 

6 9.  77 

.00931 

.060 

.096 

74.  29 

. C0911 

.080 

.131 

79.  03 

. 00869 

.105 

.17? 

83.47 

.00808 

. 130 

.216 

8 7.43 

.00742  .00208 

.00309 

.00185 

.01259 

.471 

.147 

.155 

. 764 

90.  56 

.00677  .00175 

.00202 

.00164 

.01134 

.476 

.14* 

.190 

.311 

94.55 

.00601  .00156 

.00249 

.00145 

.00986 

.482 

.147 

. 225 

.369 

97.68 

.00494  .00159 

.00219 

.00131 

.00872 

.467 

. 15  d 

.275 

. 451 

101.69 

.00391  .00120 

.00181 

.00102 

.00700 

.456 

• 1 4 < 

. 325 

.536 

104. 98 

.00306  .00097 

.00144 

.00080 

.00547 

.464 

• 1 4 E 

.375 

• MS 

1 0 7.  61 

.00236  .00082 

.00109 

.00063 

.00427 

.453 

.148 

.460 

. 738 

110.82 

.00156  .00053 

.00067 

.00040 

.00276 

• 4*0 

.14* 

.525 

.861 

113.34 

.00087  .00040 

.00043 

.00024 

.00170 

.407 

.141 

.600 

.984 

114.85 

. 00041 

ISO 


